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^ ■ Abstract 

Context. The interstellar medium is enriched primarily by matter ejected from evolved low and intermediate mass stars. The outflow 
from these stars creates a circumstellar envelope in which a rich gas-phase chemistry takes place. Complex shock-induced non- 
equilibrium chemistry takes place in the inner wind envelope, dust-gas reactions and ion-molecule reactions alter the abundances in 
I the intermediate wind zone, and the penetration of cosmic rays and ultraviolet photons dissociates the molecules in the outer wind 

• , region. 

■ Aims. Little observational information exists on the circumstellar molecular abundance stratifications of many molecules. 
Furthermore, our knowledge of oxygen-rich envelopes is not as profound as for the carbon-rich counterparts. The aim of this pa- 

JL. per is therefore to study the circumstellar chemical abundance pattern of 11 molecules and isotopologs ( "CO, CO, SiS, SiO, 

^ ! *^SiO, ^"SiO, HCN, CN, CS, SO, SO2) in the oxygen-rich evolved star IK Tau. 

Methods. We have performed an in-depth analysis of a large number of molecular emission lines excited in the circumstellar envelope 
around IK Tau. The analysis is done based on a non-local thermodynamic equilibrium (non-LTE) radiative transfer analysis, which 
calculates the temperature and velocity structure in a self-consistent way. The chemical abundance pattern is coupled to theoretical 
outer wind model predictions including photodestruction and cosmic ray ionization. Not only the integrated line intensities, but also 
the line shapes, are used as diagnostic tool to study the envelope structure. 
^ ' Results. The deduced wind acceleration is much slower than predicted from classical theories. SiO and SiS are depleted in the enve- 

' lope, possibly due to the adsorption onto dust grains. For HCN and CS a clear difference with respect to inner wind non-equilibrium 

predictions is found, either indicating uncertainties in the inner wind theoretical modeling or the possibility that HCN and CS (or 
. the radical CN) participate in the dust formation. The low signal-to-noise profiles of SO and CN prohibit an accurate abundance 

determination; the modeling of high-excitation SO2 lines is cumbersome, possibly related to line misidentifications or problems with 
the coUisional rates. The SiO isotopic ratios (^®SiO/^*SiO and ^''SiO/^^SiO) point toward an enhancement in ^*SiO compared to 
results of classical stellar evolution codes. Predictions for H2O emission lines in the spectral range of the Herschel/HIFI mission are 
performed. 



Key words. Line: profiles. Radiative transfer. Stars: AGB and post-AGB, (Stars): circumstellar matter. Stars: mass loss. Stars: 
individual: IK Tau 



I ' 1 . Introduction a C/O ratio less than unity, carbon-rich C-stars have C/O > 1, 
[ and for S-stars C/O is - 1 . 
Asymptotic Giant Branch (AGB) stars are well known to release ^ ^^^^ ^^^^^^^ ^^^^ carbon-rich stars in 
significant amounts of gas and dust in the interstellar medmm ^^-^^ ^ ^-^^ chemistry takes place. This is reflected by the de- 
via (copious) mass loss. This mass loss dommates the evolu- ^^^^-^^ ^^^^ ^-^^^^^^^ chemical compounds, including un- 
tion of the star and ultimately, when the stellar envelope is ex- ^^^^j ^^^^^^ ^^^^ ^^^-^^^^^ i n die CSE of IR C +10216, the 
hausted causes the star to evolve off the AGB into the post-AGB carbon stars (e.g. ICernicharo et al.ll2()0(ih . In con- 
phase. The outflow from these evolved stars creates an envelope, ^^^^^^ ^^^j ^q_j2 compounds have been identified in the chemi- 
which fosters gas-phase chemistry. The chemical complexity in ^^jj ^^^^ interestin g oxygen-ric h evolved stars, such as IK Tau 
circumstellar envelopes (CSEs) is thought to be dominated by ^ IZiurvs et al.l l2007h . The first observations 
the elemental carbon to oxygen ratio: oxygen-rich M-stars have carbon-bearing molecules (other tha n CO) in oxygen-rich 

AGBs were somewhat unexpected (e.g. iDeguchi & GoldsmithI 

Send offprint requests to: L. Decin, e-mail: 1985; Jewell et al. 1986). Nowadays, the formation of carbon 

Leen.Decin@sterkuleuven.ac.be molecules is thought to be the result of shock-induced non- 
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iDuari et al ] ll999h and/or a complex chemistry in the outer enve- 
lope triggered by the penetration of cosmic rays and ultra-violet 
radiation (e.g. Willacv & Millar 1997). Recently, anew interstel- 
lar molecule, PO {X ^Hr), has been detected toward the enve- 
lope o f the oxygen-rich supergiant VY CMa (iTenenbaum et al.l 
l2007h . Phosphorus monoxide is the first interstellar molecule de- 
tected that contains a P-O bond, a moiety essential in biochemi- 
cal compounds. It is also the first new species identified in an 
oxygen-rich, as opposed to a carbon-rich, circumstellar enve- 
lope. These results suggest that oxygen-rich shells may be as 
chemically diverse as their carbon counterparts. 

Circumstellar molecules have been extensively observed, 
both in the form of surveys of a single molecular species and in 
the form of searches for various molecular species in a limited 
number of carefully selected sources. The aim of these studies 
was to derive (i.) the mass-loss rate (from CO rotational lines) 
or (ii.) molecular abundances. For this latter purpose, several 
methods exist, each with varying degrees of comple xity. ( 7.) For 
example, [Bujarrabal et al. (1993) and loiofsson etall ( Il998h 
showed that simple molecular line intensity ratios, if properly 
chosen, may be used to study the chemical behaviour in CSEs. 
The use of line intensity ratios has the advantage of requiring no 
assumptions about a circumstellar model, but it also limits the 
type of conclusions that can be drawn. (2.) Several authors have 
derived new constraints on chemical and circumstellar models 
based on the simplifying assumption of unresolved optically 
thin emission thermalized at one excit ation temperature (e 
Lindqvist etal. "19881; lOmont et al.ll 19931: iBuiarrabal et"al1[T9» 
Kiin et al. 2010). (3.) Later on, observations were (re)-analyzed 
based on a non-LTE (non lo cal thermodynamic equilibrium) ra- 
diative transfer model (e g. iBieging et al.l l2000t iTevssier et al.l 
l2006HSchoier et al.l2007ah . In this study, we will go one step fur- 
ther and abandon or improve few of the assumptions still made 
in many non-LTE analyses. 



1 . Quite often, the temperature structure — being the most im- 
portant factor determining the molecular line excitation — 
is approximated usi ng a power-law (e.g. Bieging et al. 2000:; 
ITevssier et al.ll2006l) . Effects of different heating and cooling 
mechanisms are hence not properly taken into account. For 
instance, in the outermost parts of the envelope the tempera- 
ture profile deviates from a power law distribution once the 
influence of photoelectric heating by the external interstel- 
lar radiation field becomes important (e.g. Crosas & MentenI 



ll997Hjusttanont et al.lll994tlDecin et al.il2006l) 

2. The shell is often assumed to expand at a co nstant velocity 
(e. g. IBieging et alJl2000l;[Schoier et al.ll2007ah . However, for 
molecular lines primarily formed in the wind acceleration 
zone, the effect of a non-constant velocity structure on the 
derived molecular abundance may be significant. 

3. The fractional abundances are estimated to follow an expo- 
nential or Gaussian distribution, assuming that the molecules 
are formed in the inner envelope, and ph otodissociated or ab- 
sorbed onto dust grains further out (e.g. Bieging et al.ll200dt 
iGonzalez Delgado et al.ll2003l: ISchoier et al.il2007ah . The ef- 
fect of extra formation and/or depletion processes in the en- 
velope can hence not be taken into account. 

4. Often, a maximum of two molecules (CO and one other) 
is analyzed at onc e (e.g. iGonzalez Delgado et al.l 120031: 
ISchoier e t al.ll2007ah . 

5. Integrated line intensities are often used as a criterion to 
analyse the circumstellar chemical structure. However, line 
shapes provide us with strong diagnostic constraints to pin- 



point the wind acceleration, which in turn has an influence 
on the deduced fractional abundances. 

In this paper, we will study the circumstellar chemical abun- 
dance fractions of eleven different molecules and isotopologs 
in the oxygen-rich AGB star IK Tau based on the non-LTE ra - 
diative transfer code GASTRoNOoM dPecin et al.ll2006ll2007l) . 
which computes the temperature and velocity structure in the en- 
velope in a self-consistent way. Chemical abundance stratifica- 
tions are coupled to theoretical n on-equilibrium (non-TE) pre- 
dictions in the outer envelope by IWillacv & Millad (Il997h and 
co mpared to the shoc k-in duced non-TE inne r wind predictions 
bv iDuari et a D (Il999h and lCherchnei (l2006h . IK Tau has been 
chosen for study because of the wealth of observations which are 
available for this target and the fact th at its envelope is thought 
to be (roughly) spherically symmetric ( Lane et anil987t iMarvell 



120051: iHale et aUl 19971: iKim et alJl2O10l^ 

IK Tau, also known as N ML Tau, was discovered in 1965 
by iNeu^ebauer et alJ d 19651) . It is an extremely red Mira-type 
variable with spectral ty pe ranging from M8.1 to M l 1.2 and a 
period around 470 days dWing & Lockwoodlll973h . From dust 
sh ell motions de tected at 1 1 fim using the ISI interferometer. 
Ha le et al ] (119971) deduced a distance of 265 pc. This is in good 
agreement with the results of lOlofsson et al.l d 1998b who com- 
puted a distance of 250 pc from integrated visual, near-infrared 
and IRAS data using a period-luminosity relation. Estima ted 
mass-loss rates range from 3.8 x 10"^ ( Neri et al] 1 19981) to 
3 X IQ-s Mg/yr dOonzalez Delgado et al.ll2003l) . IK Tau's prox- 
imity and relatively high mass-loss rate (for a Mira) facilitates 
the observation of molecular emission lines. 

In Sect. |2] we present the molecular line observational data 
used in this paper. Sect. [3] describes the background of the ex- 
citation analysis: the radiative transfer model used, the molec- 
ular line data and the theoretical ideas on molecular abundance 
stratification in the envelope. Sect. |4] describes the results: we 
first focus on the velocity structure in the envelope with special 
attention to the acceleration zone, after which the derived stel- 
lar parameters are discussed. Thereafter, the abundance structure 
for each molecule is derived and compared to the theoretical in- 
ner and outer wind predictions and observational results found 
in the literature. The time variability and SiO isotopic ratios are 
discussed in Sect. |5] and water line predictions are performed in 
Sect. 15.31 We end with some conclusions in Sect.|6l 



2. Observational data 

Part of the observations were obtained from our own observ- 
ing programs scheduled at the JCMT, APEX and IRAM. These 
observations and the data reduction are described in Sect. 12.1! 
Other data are extracted from the literature and summarized 
briefly in Sect. 12.2! An overview is given in Table [1] 

2. 1. Observations and data reduction 

The i2cO(2-l), i2cO(3-2), i2cO(4-3) and the i3cO(2-l) ob- 
servations were extracted from the jcmtQ archive. Additional 
data with the APEXEI 12 m telescope were obtained for the 
i2C0(3-2), i2cO(4-3), i2cO(7-6), and i3cO(3-2) molecular 



The James Clerk Maxwell Telescope (JCMT) is operated by The 
Joint Astronomy Centre on behalf of the Science and Technology 
Facilities Council of the United Kingdom, the Netherlands Organisation 
for Scientific Research, and the National Research Council of Canada 

~ APEX, the Atacama Pathfinder Experiment, is a collaboration 
between the Max-Planck-Institut fur Radioastronomie, the European 



2 



L. Decin et al.: The circumstellar chemistry in the O-rich AGB IK Tau 



line transitions. With the 30 m telescope of the Institut de Radio 
Astronomie Millimetrique (IRAMj^ molecular line observa- 
tions were performed in December 2006. During this observing 
run, data on the CO(2-l), SiS(8-7), SiS(12-ll), SOzCWa n - 

142,12), S02(43,l - 42,2), S02(33,l " 82,2), S02(53,3 " 52,4), 

and HCN(3-2j3 line transitions were obtained. 

The JCMT and APEX observations were carried out using 
a position-switching mode. The IRAM observations were done 
in the wobbler-switching mode with a throw of 60". The fre- 
quency resolution for the JCMT-data equals 0.0305 MHz, for the 
APEX data it is 0.1221 MHz. The resolution was 1.25 MHz for 
the 3 and 2mm IRAM observations and 1 or/and 4 MHz at 1.3 
and 1 mm, resulting in resolution slightly higher than 1 km/s for 
these observations. 

The JCMT data reduction was performed with the SPLAT 
devoted routines of STARLINK, the APEX and IRAM data were 
reduced with CLASS. A polynomial of first order was fitted to an 
emission free region of the spectral baseline and subtracted. To 
increase the signal-to-noise ratio, the data were rebinned to a res- 
olution of ^ 1 km/s so that we have at least 40 independent res- 
olution elements per line profile. The antenna temperature, TJ, 
was converted to the main-beam temperature (Tmb = r^/^ymb), 
using a main-beam efficiency, 77,„b as specified in Table [T] The 
absolute uncertainties are ^20 %. 



2.2. Literature data 

To have better constraints on the chemical abundance pattern in 
the wind region around IK Tau, additional data were taken from 
the literature (see Table [T]i. 

High-quality observations were performed by iKimI (|2006|) 
with the APEX telescope in Chile during observing periods 
in November 2 005, April 2006 and August 2006 (see also 
lKimetalJ2010l hereafter referred to as Paper I). In total, 34 tran- 
sitions from 12 molecular species, including a few maser lines, 
were detected toward IK Tau. 

ISchoier et al.l ( l2007bl) published the observations of four SiS 
lines: the (5-4) and (6-5) rotational line transitions were ob- 
tained with the Onsala Space Observatory (OSO) telescope, the 
(12-11) and (19-18) rotational line observations were performed 
with the JCMT telescope. 

SiO thermal radio line emission from a large sample 
of M-type AGB stars, including IK Tau, was studied by 
iGonzalez Delgado et all (l2003h . The SiO(2-l) line transition 
was obtained with the OSO telescope, the SiO(5^) and (6-5) 
transitions with the Swedish ESO Submillimeterwave Telescope 
(SEST). 

^^CO line data were obtained bv lTevssier et alj ( l2006l) with 
th e IRAM and JCMT tele scope. The CO(l-O) fine was observed 
bv lRamstedt e t al.' (2008) in December 2003 with the 20 m OSO 
telescope, loiofsson et aLa998b reported on the detection of the 
CS(2-1) fine at 98 GHz with the OSO telescope with an inte- 
grated intensity of 0.5 K km/s. 



Table 1. Overview of the molecular line transitions used in 
this research, with indication of the frequency, the upper en- 
ergy level, the telescope, the main beam half power beam width 
(HPBW) and main beam efficienciy (r/,„b)- 



Southern Observatory, and the Onsala Space Observatory. Program IDs 
are077.D-0781 and 077.D-4004. 

^ IRAM is supported by INSU/CNRS (France), MPG (Germany), 
and IGN (Spain). 

When the isotopic notation is not given, the molecular line transi- 
tion is from the main isotopolog. 
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3. Excitation anaiysis 

3.1. Radiative transfer model 

The observed molecular line transitions provide information 
on the thermodynamic and chemical structure in the enve- 
lope around IK Tau. The line profiles were modeled with 
our non-LTE (non-Local Thermo dynamic Equilibri um) radia- 
tive transfer code GASTRoNOoM ( iDecin et alj|2006h . The code 
(1) calculates the kinetic temperature and velocity structure 
in the shell by solving the equations of motion of gas and 
dust and the energy balance simultaneously; then (2) solves 
the radiative transfer equation in the co-moving frame using 
the Approximate Newton-R aphson operator as developed by 
ISchonberp & HernpS (Il986h and computes the non-LTE level- 
populations; and finally (3) determines the observable line pro- 
file by ray-tracing . For a full description of the code we refer to 
IDecin et al.l(l2006h . 

The main assumption of the code is a spherically symmetric 
wind. The mass-loss rate is allowed to vary with radial distance 
from the star. The local line width is assumed to be described by 
a Gaussian and is made up of a microturbulent component with 
a Doppler width of 1.5km/s and a thermal component which is 
calculated from the derived kinetic temperature structure. 

Two major updates have been made since the original publi- 
cation in IDecin etall ( l2006l) . 

- The code now iterates on steps (1) and (2) to obtain the ki- 
netic temperature structure in a self-consistent manner from 
solving the energy balance equation, where the CO and H2O 
line cooling (or heating) are directly obtained from the exci- 
tation analysis, i.e. 



(1) 



I U>1 



where ni and n„ are the level populations in the lower and 
upper levels participating in the transition at rest frequency 
i^ui, and Jul and jiu are the CO-H2 collisional rate coeffi- 
cients. The cooling rate A (in erg s^^ cm^""*) is defined as 
positive for net cooling. In case of IK Tau the water line cool- 
ing dominates the CO line cooling by more than an order 
of magnitude in the inner wind region; for regions beyond 
10^^ cm the CO line cooling dominates over H2O cooling 
with the adiabatic cooling being the dominant coolant agent. 
While the original version of the code approximates the stel- 
lar atmosphere with a blackbody at the stellar effective tem- 
perature, an additional option is now implemented to use 
MA RCS theoretical model a t mospheres and theoretica l spec- 
tra dGustafsson et al.l 120081: iDecin & Erikssonll2007l) to es- 
timate the stellar flux. Molecular species in the CSE, less 
abundant than CO and with larger dipole moments are pri- 
marily excited by infrared radiation from the c entral star 
(with the possible exception of HCN, lJuralll983l and H2O, 
see Sect. 15.3b . For CO, the infrared radiation competes with 
rotational excitation by collisions and by trapped rotational 
line photons to determine the populations of the rotational 
levels (Knapp & Morris 1985). For those minor species, the 
blackbody approximation of the stellar flux may lead to in- 
accurate absolute intensity predictions in the order of 5 to 
20%. 



a proper determination of the envelope structure. This is par- 
ticularly true since (i.) the line shapes contain valuable infor- 
mation on the line forming region, e.g. a Gaussian line profile 
points toward line formation partially in the inner wind where 
the stellar winds has not ye t reached its full terminal velocity 
(iBuiarrabal & Alcolealll99ll) . and (ii.) while the absolute uncer- 
tainty of the intensity ranges between 20 to 50 %, the relative 
accuracy (of the line shapes) is of the order of a few percent. 
In contrast to most other studies, we therefore will not use a 
classical x^-analysis using the integrated line intensity, but will 
perform the model selection and the assessment of goodness-of- 
fit using the log-likelihood fun ction on the full line profiles as 
described in lDecin etaD(l2007h . 



3.2. Molecular line data 

In fliis paper, line transitions of CO, SiO, SiS, CS, CN, HCN, SO, 
and SO2 will be modeled, and H2O line profile predictions for 
the Herschel/HIFI instrument will be performed. The molecular 
line data used in this paper are described in AppendixlAl 



3.3. Molecular abundance stratification 

Theoretical chemical calculations clearly show that the frac- 
tional abundances (relative to H2) vary throughout the envelope. 
Chemical processes responsible for the molecular content are 
dependent on the position in the envelope (see Fig. [T]i. In the 
stellar photosphere and at the inner boundary of the envelope, 
the high gas density and temperature ensure thermal equilibrium 
(TE). TE is suppressed very close to the photosphere because 
of the action of pulsation-driven shocks propagating outwards. 
Furthermore the regions of strong shock activity correspond to 
the locus of grain formation and wind acceleration. This region 
is referred to as the inner envelope (or inner wind) which extends 
over a few stellar radii. At larger radii (-^5 to 100 R*) the newly 
formed dust grains interact with the cooler gas. Depletion or for- 
mation of certain molecular/atomic species may result from such 
an interaction and these layers are referred to as the intermedi- 
ate envelope. This is also the region where parent molecules, 
injected in the envelope, may begin to break down, and daugh- 
ter molecules are formed. At still larger radii (>100R*), the 
so-called outer envelope is penetrated by ultraviolet interstellar 
photons and cosmic rays resulting in a chemistry governed by 
photochemical processes. 

Since our modeling results will be com pared to chemical 
abundance predictions in the outer envel ope by Willacv & Millarl 
(1 19971) and in the inner envelope by iDuari et al.l d 19991) and 
ICherchnefS (l2006h . we first briefly describe these studies in 
Sect. l3.3.Tl and Sect. |332] respectively. In Sect. 13.41 we discuss 
how we have implemented this knowledge in the modeling of 
the molecular line transitions. 



IDecin et al.l (l2006i l2007l) and ICrosas & MentenI (Il997h 
demonstrated that not only the integrated line intensities, but 
particularly the line shapes should be taken into account for 



3.3.1 . Chemical stratification in the outer envelope 

Th e chemistry in the outer e nvelope of IK Tau has been modeled 
bv lWillacv & Millaij (Il997h . This chemical kinetic model aims 
at deriving the abundance stratification in the outer envelope (be- 
tween 2 X 10^^ and 2 x 10^® cm). The chemistry is driven by a 
combination of cosmic -ray ionization and ultraviolet radiation 
and starts from nine parent species injected into the envelope 
(see Table|2]i. The CSE of IK Tau was assumed to be spherically 
symmetric with a constant mass loss rate and a constant expan- 
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Figure 1. Schematic overview (not to scale) of the circumstellar envelope (CSE) of an oxygen-rich AGB star. Several chemical pro- 
cesses are indicated at the typical temperature and radial distance from the star in the envelope where they occur. The nomenclature 
as used in this paper is given. 



sion velocity of 19km/s. The temperature was described by a 
power law 

r(r)=T(ro)(^)"' (2) 

with ro = 2 X 10^^ cm and T(ro) taken to be 100 K for M> 
5 X IQ-'^Mo/yr and 300 K otherwise. For IK Tau T{ro) was 
assumed to be 100 K. The derived fractional abundances for the 
molecules studied in this paper are represented in Fig.|2] 

Table 2. The fra ctional abundance ( relativ e to H2) taken for the 
parent species bv lWillacv & Millad (119971) . 



Parent 
species 


Fractional 
abundance 




Parent 
species 


Fractional 
abundance 




CO 


4.00 X IQ- 


4 


N2 


5.00 X 10" 


5 


H2O 


3.00 X 10" 


4 


H2S 


1.27 X 10" 


5 


CH4 


3.00 X IQ- 


5 


SiS 


3.50 X 10" 


6 


NH3 


1.00 X IQ- 


6 


SiO 


3.15 X 10" 


5 


He 


1.00 X 10" 


1 


PH3 


3.00 X 10" 


8 


HCl 


4.40 X 10" 


7 









The model of IWillacv & Millaj (Il997h succeeded in repro- 
ducing the observed values of certain species but failed for 
some other molecular abundances: the calculated abundance of 
HCN was too low and the injected abundance o f the parent 
specie s SiS was about 10 times higher than observed. iDuari et"al1 
(Il999h noted that the input molecular abundances of some par- 
ent species are sometimes questionably high since there exists no 
observational or theoretical evidence for the formation of these 
species in the inner and intermediat e envelopes of O-ri ch Miras 
(see Sect. I3.3.2l i. More importantlv. iDuari et al. l ([T999l) showed 
that HCN should form in the inner envelope or extended stellar 
atmosphere due to non-equilibrium shock chemistry, and may be 
a parent species injected to the outer envelope. Recent observa- 
tional stu dies also indicate that HCN must be formed in the inner 
envelope dBieging et al. 2000; Marvel 2005). These resuhs are 
in contrast to the modeling efforts of Willacy & Millaj (Il997h 
where HCN was not yet considered as a parent species. 



3.3.2. Chemical stratification in the inner envelope 

Carbon-bearing molecules have been identified in the envel opes 
of many oxygen-rich AGB stars (e.g. jBuiarrabal et al.ll994l) and 
it was first thought that the observed carbon species were pro- 
duced in the out er wind of O-rich s tars via photochemical pro- 
cesses. However. iDuari et"al ] (ll999h showed that shock-induced 
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Figure 2. Theoretically predicted molecular abundance stratifi- 
cation in the envelo pe of IK Tau. The full l ines represent the 
results as derived by I Willacy & MiUd ([T997h for the outer en- 
yelope of IK Tau assuming a mass-loss rate of 1 x 1 0~^ M( T,/yT. 
The di amonds represent the abundances deriyed by Cherchneff 
(l2006h in the framework of non-equilibrium chemistry due to the 
passage of shocks in the inner eny elope of TX Cam a t 2 i?*; the 
filled circles indicate the results by lDuari et al ] (ll999h forIKTau 
at 2.2 Ri,. The inner wind predictions are projected onto the dust- 
condensation radius deriyed in Sect. 14.11 The dotted line repre- 
sents the Gauss i an dis tribution of circumstellar H2O as used by 
lMaerckeretal1 ( l2008l) . 



non-equilibrium chemistry models predict the formation of large 
amounts of a few carbon species, such as HCN, CS and CO2, in 
the inner enyelope of IK Tau: these molecules are hence formed 
in the post-shocked layers and are then ejected in the outer wind 
as 'parent' species. For some parent species, the non-equilibrium 
chemistry does not alter significantly the initial photospheric TE 
abundances. But other species, abundant in the TE photosphere, 
are quickly destroyed in the outflow by the non-equiblibrium 
chemistry generated by shocks (e.g., OH, SiS and HS). Again 
other species (such as SO) appear to be absent in the inner 
regions of the wind, and are thought to be produced by ion- 
molecule reactions in the photo-dissociation regions of the outer 
wind. 

ICherchnel(l2006h continued on the study of shock-induced 
non-equilibrium chemistry in the inner wind of AGB stars. She 
demonstrated that whateyer the enrichment in carbon of the star 
(i.e. the C/O ratio), the atomic and molecular content after the 
passage of the first shock in the gas layers just aboye the stellar 
photosphere is yery much the same, and in many cases, totally 
different from what would be expected from thermodynamic 
equilibrium (TE) calculations. In the case of the oxygen-rich en- 
yelope around TX Cam — being almo st a stellar twin of I K Tau, 
but with slightly lower luminosity — ICherchnefB (l2006l) found 
that while, e.g., the TE abundance of HCN (CS) is predicted to 
be ~1.9 X 10"" (-2.5 x 10""), the non-TE fractional abun- 
dances at 2.5 R* are predicted to be ~9 x 10~^ (~1 85 x 1 0~^). 
The fractional abundances deriyed by I Cherchneff (l2006h dif- 
fer from the abundances of the injected parent molecules in 
the study of iWill a cy & M illar (1997) (see Table |2] and Fig. |2]i: 
IWillacy & MiHarl (Il997k did not consider CS and HCN as be- 
ing parent molecules, and the (injected) abundance of SiS (3.5 x 



10 ^) is much highe r than the abundance stratification deriyed 
by lCherchnei(l2006h (see their Fig. 8). 

3.4. Modeling strategy 

3.4.1 . Envelope structure as traced by the CO lines 

The physical properties of the circumstellar gas, such as the 
temperature, yelocity and density structure, are determined 
from the radiatiye transfer modeling of the multi-transitional 
(sub)millimetre CO line obseryations. Since higher- J lines are 
formed at higher temperature, different transitions offer the pos- 
sibility to trace different regions in the enyelope. The highest 
CO energy leyel traced is the CO J"p = 7 leyel at 154.8 K. The 
ayailable rotational CO lines will hence be good tracers for the 
region beyond ~ 1 00 R^, , but they do not put strong constraints on 
the temperature in the inner CSE. The upcoming Herschel/HIFI 
mission will be crucial in the study of the temperature structure 
in this inner wind region. 

An extensiye grid has been calculated with parameters rang- 
ing from 2000 to 3000 K for the stellar temperature Tcs, from 
1 X lO^'' to 6 X lO^'^ cm for the stellar radius R^,, an inner (dust 
condensation) radius between 2 and 30i?,t, distance between 
200 and 300 pc, and a constant mass-loss rate between 1 x 10^^ 
and 5 x 10"^ Mpi lyr. As briefly ex plained in Sect. 13.11 a log- 
likelihood method (iDecin et al.l l2007!) is used to find the best-fit 
model and deriye a 95 % confidence interyal for the model pa- 
rameters. The results will be presented in Sect.|4] 

3.4.2. Abundance stratification through the envelope 

From the descriptions of theoretical abundance estimates in the 
inner and outer enyelope (Sect. l33n - I3.3.2l i it is clear that there 
is still some debate about the abundance structure in the enye- 
lope. SiS and HCN were already giyen as an example, but other 
molecules such as e.g. CS and SO also pose a problem. This 
is illustrated in Fig. |2] where one notices for a few molecules 
a large difference between the theoretically predicted fractional 
abundance i n the inner enyelope by Duari et al. (199^ and 
ICherchneS (l2006l) and the abun dance of the parent mole cules 
injected in the outer enyelope by I Willacy & MiliaJ (Il997l) . One 
of the big questions still existing concerns the modifications of 
the molecular abundances in the intermediate wind region due to 
gas-grain reactions. Currently, no theoretical efforts haye been 
made to model this region in terms of molecular 'leftoyers' after 
the dust formation has occurred. In the case of O-rich enyelopes, 
it is thought that CO, HCN and CS are quite stable and trayel the 
entire enyelope unaltered until they reach the photo-dissocation 
region of the outer wind as these molecules do not participate 
in the formation o f dust grains such as silicates and corundum 
(iDuari et al. l [T999h . In contrast, SiO is a candidate molecule for 
depletion in the intermediate wind region due to the formation 
of Si02 (yia a reaction with OH) whose condens ation product, 
silica, is tentatiyely identified in post-AGB stars jMolster et al.l 
I2002h and is claimed to b e the carrier of the 1 3 ym feature in low 
mass-loss rate AGB stars dSpeck et al.l2000l but other studies ar- 
gue that this feature is d ue to spinel) . The th eoretical modeling of 
iDuari et all d 19991) and ICherchnel (l2006l) predict an SiS abun- 
dance 2 to 3 orders of magnitude lower than the obseryed yalue, 
indicating that SiS is produced in the o uter enyelope of IK T au. 
Howeyer, recent obseryational results by lDecin et al ] (l2008bh ar- 
gue for a formation process in the inner enyelope. 

From the aboye arguments, it is clear that we should allow 
some yariation in modeling the abundance structure in the enye- 
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lope. However, one should also realize that we sometimes only 
have 2 rotational transitions of one isotopolog at our disposal 
with a restricted range in excitation temperature. The highest up- 
per level energy traced is the SiS(20-19) transition at ~^183K; 
hence none of the studied transitions is sensitive to the abun- 
dance in the inner envelope (R^ 5 i?^). In order to use some 
prior knowledge on the (theoretical) photo-dissociation rate in 
the outer regions and to allow for a depletion or an extra for- 
mation process in the intermediate/outer envelope, we therefore 
have opted to divide the envelope in different regimes (see also 
Tableland FiglTJ: (i.) in the dust-free zone {R < i?innor) the 
abundance is constant (/i(mol)), (ii.) between i?inncr and i?niax 
the abundance can decrease/increase from /i(mol)(i?i,inor) to 
/2(mol)(i?inax) linearly on a log-log scale, where both i?,nax 
and /2(mol) are free parameters, (Hi.) from i?inax onwards, the 
ab undance stratificat i on fo llows the (photodissociation) results 
of lWillacv &Millaj (Il997h scaled to /2(mol) at i^^ax- In that 
way, 3 parameters (/i(mol), /2(mol), and i?max) have to be es- 
timated to determine the abundance stratification of a species. 



Table 3. Modeling assumptions of the abundance stratification. 



Region 


Variable 


Comment 


R ^ -dinner 


/o(mol) 


not sensitive 


R — -dinner 


/i(mol) 


/o(mol)= /i(mol) 




/2(mol) 


allow for depletion (e.g. due to dust 




R2 


formation) or extra formation process 


dinner < R < R2 




linear interpolation in log(R)-log(f) 


R> R2 




theor. pred. of Willacy scaled to /2(mol) 



Most studies use the photodissociation results of 
iMamon et al.l (Il988h to describe the CO spatial variation 
in the outer envelope. For other molecules, the abundance 
pattern is often assumed to be described by a simple 



Gaus si an or expontential di s tributi o n (e.g. Bieging et alj 
|200(A iGonzalez Delgado et al.l l2003t ISchoier et alJ l2007bh . 
The e-folding radius then describes the photodissociation by 
ambient UV photons penetrating the dusty envelope or depletion 
of the molecules from the gas into dust grains in the outflowing 
stellar wind. That way, however, all molecules are assumed to 
be created in the extended atmosphere or inner wind region. 
Moreover, a combination of depletion and photodissociation or 
extra depletion/formation processes in the intermediate/outer 
region can not be captur ed, and one can not use the results by 
IWillacv & MiiiaJ (1 1 9971) describing extra formation of a few 
molecules by ion-ion reactions in the outer wind region. The 
methodology outlined above (Table O captures these flaws, and 
may serve to considerably strengthen our knowledge on the 
abundance stratification in the envelope. 

As was already alluded to in the previous paragraph, the 
line profiles in this study are not sensitive to a change in abun- 
dance in the inner wind region (i? < i?i„nor)- To assess the 
abundance stratification in this reg ion, on e either needs high- 
resolution near-infrared (see, e.g., iDecin et al. 2008a) or far- 
infrared spectroscopy (as will be provided by the Herschel/HIFI 
instrument). Nonetheless, the derived abundance stratifications 
will be compared to th e theoretica l inner wind predictions by 
iDuari et al ] (Il999h and lCherchnel(l2006h . since this compaii- 
son may yield hints on the (un)reactivity of the molecules in the 
dust forming region and on uncertainties in the inner wind pre- 
dictions. 



4. Results 

Using the log-likelihood method, the parameters for the model 
yielding the best-fit to the CO line profiles are derived (see pa- 
rameters listed in the 2nd column in Table [H 'model 1'). The 
CO lines, however, only trace the envelope beyond '^lOOi?*. 
One therefore should use other molecules to put constraints on 
the structure in the inner wind region. HCN is the only molecule 
for which we have observational evidence that it is formed in 
the inner wind region: using interferometric data Marvel (2005|) 
deduced a maximum size for the HCN distribution of 3.85"(in 
diameter), or a radius of 7.2 x 10^^ cm at 250pc. They con- 
cluded that the deduced size indicates a shock origin for HCN 
close to the star and a radius for the HCN distribution lim- 
ited by photodissociation. The HCN line profiles (Fig. |7]l are 
clearly Gaussian indicating a line formation (at least partly) in 
the inner wind region, where the wind has not yet reached its 
terminal velocity. That way, HCN observations yield important 
clues on the velocity structure in the inner wind region. Using 
the stellar parameters as given in the 2nd column in Table |4] 
('model 1'), we were unable to derive a HCN-abundance struc- 
ture yielding a satisfactory fit to the line profiles. While the inte- 
grated intensities could be well predicted, the line profiles were 
flat-topped with a FWHM (full width at half maximum) being 
too large. The only way to reconcile this problem was (1) con- 
centrating the HCN abundance in the inner 2 x 10^^ cm with 
[HCN/H] = 9 X 10~^, or (2) aUowing for a smoother velocity 
law. While th e abundance in the former solution is within the 
predictions of ICherchneffI (l2006h , the ang ular distance is much 
smaller than the 3.85"diameter deduced bv lMarvell(l2005h . 



4.1. Velocity Structure 

The expansion velocity of SiO, H2O, and OH masers can be used 
to put further constraints on the velocity structure, and specifi- 
cally on the acceleration in the inner wind region (see Fig. O. 
It is clear that the velocity structure as derived from the param- 
eters of the best-fit model only based on CO fines (model 1), is 
far too steep in the inner wind region as compared to the velocity 
indications of the maser lines. This problem can be solved by ei- 
ther increasing the dust condensation radius or by allowing for a 
smoother veloci ty profile. This latter can be simulated using the 
classical /3-law (iLamers & Cassineliilll999l) with /3 > 0.5 (see 
Fig.EJ 



v{r) ~vo + {voo - Wo) ( 1 - ^ 



(3) 



with f the velocity at the dust condensation radius. 

One should realize that several assumptions are inherent to 
the velocity structure derived from solving the momentum equa- 
tion: (i.) All dust species at all different grain sizes are as- 
sumed to be directly formed at the dust co ndensation radius 
dinner- Howcvcr, theoretical results from, e.g. jGail & Sedknavj 



(Il999h show that formation and growth of (silicate) dust grains 
typically occur between 1100 and 900 K, i.e. extending over 
a few stellar radii, (ii.) The extinction efficiencies used in the 
GASTRoNOoM-code represent the Fe-rich silicate MgFeSi04. 
Thanks to their high absorption efficiencies at optical and near- 
infrared wavelength Fe-rich sili cates hke MgF eSi04 (and solid 
Fe) are efficient wind drivers (IWoiticdl2006h . However, other 
oxides or pure silicates such as AI2O3, Si02, Mg2Si04 and 
MgSiOa have low absorption efficiencies at optical and near- 
infrared wavelengths, resulting in a negligible radiative pres- 
sure on all glassy condensates. If these latter molecules were the 
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Figure 3. Velocity profile of IK Tau. Velocity data are obtained 
from mappin g of maser emis sion: SiO (Boboltz & Diamondl 
I2005h . H2O feains et al.ll20()3h . and OH (iBowers et al.lll989l) . 
The CO expansion velocity derived from our CO data is also 
indicated. The expansion velocity deduced from the CO data 
alone (see 'model 1' in Table|4|i is plotted as full black line. The 
dotted black line indicates the velocity structure taking a tur- 
bulent velocity (of 1.5km/s) into account. The green line gives 
the expansion (+turbulent) velocity deduced from both the CO 
and HCN lines (model 2 in Table |4]i. The dashed blue line rep- 
resents an even smoother expansion (H-turbulent) velocity struc- 
ture, applying Eq.[3] with /3=l (model 3 in Table|4]i. The vertical 
dashed black and green lines indicate the dust condensation ra- 
dius Rinncr- Notc that /3=0.5 in the inner wind region. 



Using the same stellar parameters as in 'model 2', but simu- 
lating a smoother velocity law still compatible with the maser 
line mapping (/3 = 1) decreases the mass-loss rate to 8 x 
10~^ AfQ/yr (model 3 in TablelH, due to the condition of mass 
conservation (model 3 in Table|4|. The narrow Gaussian line pro- 
files of the HCN lines (Sect. l43!2] i give more support to model 3 
than to model 2, for which reason the thermodynamic structure 
as deduced from model 3 (Fig.|4]i will be used to model the other 
molecular line transitions. 



4.2. Stellar parameters derived from the CO and HCN 
lines 

The stellar parameters for the best-fit model (model 3) are listed 
in Table H) The outer radius of CO is computed using the re- 
sults of Mamon et al ] (ll988h . The derived 95 % confidence in- 
tervals in Table |4] are statistical uncertainties, which should be 
interpreted in the light of the mode l assumption s of a s pherically 
symmetric wind. As described in iDecin et aP (l2007h . the log- 
likelihood function can also be used to compare different mod- 
els with a different number of parameters. In case of IK Tau, we 
have assessed the likelihood preference of a model with constant 
mass-loss rate compared to a model with mass-loss rate varia- 
tions. The preferences pointed towards the simpler model, i.e. 
with a constant mass-loss rate (of 8 x 10^^ AIq/yi). We also de- 
rived the dust-to-gas mass ratio from the amount of dust needed 
to drive a wind at a terminal velocity of 17.7 km/s for a gas mass- 
loss rate of 8 x 10~^ MQ/yr with the deduced velocity profile. A 
dust-to-gas mass ratio of 1.9 x 10^^ (or a dust mass-loss rate of 
1.52 X 10~^ Mq/yt) is obtained for model 3, with an estimated 
uncertainty of a factor ^5. 



most abundant in the envelope of IK Tau, the wind acceleration 
would be much lower No medium resolution infrared (from the 
Infrared Space Observatory - Short Wavelength Spectrometer or 
the Spitzer - InfraRed Spectrograph) data are, however, avail- 
able for IK Tau, hence we were unable to study the circum- 
stellar dust composition in detail. (Hi.) 'Complete momentum 
coupling' is assumed. This means that the grain motion every- 
where in the flow can be computed by equating the local radia- 
tive and collisional drag forces implying that virtually all of the 
momentum gained by a grain through the absorption of radia- 
tion from the stellar photosphere is transfen'ed vi a collisions to 
the atmospheric gas jMacGregor & Stencellll992l) . For physical 
conditions typical of the circumstellar en velopes of oxygen-rich 
red giants. iMacGregor & Stencel ( Il992h find that silicate grains 
with initial radii smaller than about 5 x 10^^ fim decouple from 
the ambient gas near the base of the outflow, (iv.) The momen- 
tum equatio n used in the GASTRoNOoM-code (see Eq. 2 in 
iDecin et al]|200 6). implicitly assumes that the mass outflow is 
steady in time and tha t the circumstellar dust is o ptically thin 
to the stellar radi ation (iGoldreich & Scovillelll976l) . Dust emis- 
sion modeling bv lRamstedt et 2008 ) suggests a circumstellar 
envelope which is slightly optically thick at 10 fim (tiq = 1.2). 
These results suggest that the acceleration of the gaseous parti- 
cles in the inner wind might be slower than deduced from solving 
the momentum equation ('model 1' and 'model 2'), since not all 
dust species take part in the momentum transfer. 

Solving the momentum equation and taking both the CO and 
HCN line profiles into account, the inner radius of the dusty en- 
velope is shifted towards higher values, i?iii„cr = 1.3 X 10^"* cm, 
and the mass-loss rate slightly increases ('model 2' in Table|4|. 



4.3. Fractional abundances 

Using the thermodynamic envelope structure derived above (see 
Fig. IHi, the abundance stratification of all molecules is derived. 
A comparison to the theoretical inner and outer wind predictions 
(as discussed in Sect. l3.3.Tl43.3.2b is given in Fig.|5] The studied 
molecular line transitions are not sensitive to the full envelope 
size, but have a limited formation region. The part in the enve- 
lope we can trace by combining the different available rotational 
line transitions is indicated with vertical dashes in Fig. |5] and 
tabulated in Table |2l 

A comparison between observed and predicted line profiles 
and a discussion on the deduced abundance stratification are 
given for each molecule separately in the following subsections. 
We will always first briefly describe the deduced abundances, 
then compare the results to the theoretical inner and outer wind 
predictions and finally compare to other results found in the liter- 
ature (see Table|6]l. For the literature results, a difference is made 
between studies based on the assumptions of optically thin unre- 
solved emission and a population distribution thermalized at an 
excitation temperature that is constant throughout the envelope, 
and those based on a full non-LTE radiative transfer calculation. 
One should also realize that most studies make use of integrated 
line intensities, and do not deal wi th a full line profile a nalysis as 
is done here. With the exception of lOmont et al.l(ll993h . the other 
literature studies listed in Table |6] assume the shell to expand at 
a constant velocity. As discussed in Sect. 14.11 the Gaussian line 
profiles of the HCN and a few of the SiO lines are the result 



Estimates of these ranges are found by considering the place where 
I{p)p^, with / the intensity and p the impact parameter, is at half of its 
maximum value. 
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Table 4. Parameters of the models with best goodness-of-fit for IK Tau. Details for each model are given in the text. The numbers 
in italics are input parameters that have been kept fixed at the given value. 





model 1 


model 2 


model 3 


95% confidence 


based on 


CO 


CO & HCN 


CO & HCN 


interval for 


v{r) 


momentum eq. 


momentum eq. 


power-law 13 — 1 


model 3 


parameter 


value 


value 


value 





Teff[K] 


2200 


2200 


2200 


100 


7?, [10" cm] 


2.5 


1.5 


1.5 


1.0 


[CO/H2] [10-*] 


2 


2 


2 




distance [pc] 


250 


265 


265 


20 


-dinner [^*] 


2.2 


8.7 


8.7 


2 


-Router [R*] 


15170 


27000 


27000 




Voo [kms"^] 


17.7 


17.7 


17.7 




M{r) [Mg/yr] 


8 X 10-*^ 


9 X lO"** 


8 X lO"*^ 


1 X 10 


i^C/^C 


20 


16 


14 


2 





Figure 4. Thermodynamic structure in the envelope of IK Tau as derived from the ^^CO and HCN rotational line transitions for the 
stellar parameters of model 3 in Table |4] Upper left: Estimated temperature profile, upper right: estimated gas and drift velocity, 
lower left: cooling rates, and lower right: heating rates. The start of t he dusty envelope, Rjnne , is indicated by the dotted line. The 
temperature structure as used in the chemical kinetic calculations by IWillacv & Millar! (Il997h is plotted in the upper left panel in 
dashed lines. 



of line formation partially in the inner wind, where the stellar vational studies assuming a constant expansion velocity, will be 
wind has not yet reached its full terminal velocity. As a result, unable to predict the line profiles properly, 
more emission is produced at velocities near the line center than 
would be the case for a uniform-velocity wind. Hence, obser- 
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C-bearing molecules Si and S bearing molecules 




1 10 100 1000 10000 1 10 100 1000 10000 

radius [R.] radius [R.] 



Figures. Predicted abundance stratifications [mol/Htot] (full line, see Fig.|2]l are compared to the deduced abundance structures 
(dashed lines). In the left panel, the carbon bearing molecules are shown, the right panel gives the Si and S bearing molecules. For 
CN, an alternative solution is given in dotted lines. For each molecule, the line formation region traced by the observed molecular 
lines is indicated by the vertical dashes (see also Table|5]). 



Table 5. Molecular fractional abundance relative to 
Htot=n(H)+2n(H2) (seeFig.|5j. 





/i(mol) 




/2(mol) 




R2 [R.] 


ranges [i?*] 




1.0 X 10 


-4 


1.0x10" 


-4 


25170 


100-25000 




7.1x10" 


6 


7.1x10" 


^6 


25170 


500-15000 


Sis 


5.5x10" 


6 


4.0x10" 


-9 


120 


40-80, 1500-6000 


2«SiO 


8.0x10" 


6 


2.0x10" 


-7 


180 


70-10000 


29SiO 


3.0x10" 


7 


7.5x10" 


-9 


180 


70-2000 


30SiO 


1.0x10" 


7 


2.5x10" 


-9 


180 


70-2000 


cs 


4.0x10" 


8 


4.0x10" 


-8 


1160 


300-4000 


so 


2.0x10" 


7 


1.0x10" 


-9 


400 


200-3000 


S02 


1.0x10" 


6 


1.0x10" 


-7 


1000 


50-400 


HCN 


2.2x10" 


7 


4.0x10" 


-9 


500 


50-500 


CN 


3.0x10" 


8 


3.0x10" 


-8 


1000 


500-3500 








alternative solution 




CN 


1.0x10" 


10 


1.0x10" 


-10 


450 


900-3500 



4.3.1. CO 

Results: A comparison between the observed rotational ^^CO 
and ^■^CO lines and theoretical predictions is shown in Fig. |6] 
The ^^CO and ^■^CO lines are very well reproduced by the 
GASTRoNOoM-predictions, both in integrated intensities and 
in line shapes. Only the IRAM i2cO(l-0) and i2cO(2-l) Hnes 
are slightly over-predicted. It is, however, not the first time that 
the non-c ompatibility of the IRAM absolute flux level is reported 
(see, e.g. lDecin et al.ll2008aHSchoier et al.ll2006l) . 



Comparison to theoretical predictions: The CO fractional 
abundance assumed in all observational studies (see also 
Table |6]l is always lo wer than th e dedu ced i nner wind theoreti - 
cal non-TE values of iDuari et al ] (119991) and lCherchneffl(l2006h . 
The non-TE theoretic al values are comp arable to the TE- value of 
6.95 X 10"-* at 1 (iDuariet al.lll999t) . Increasing the CO frac- 
tional abundance by a factor 5, would decrease the mass-loss 
rate by a factor ~2.8 to reproduce the observed CO rotational 
line p rofiles (using the scaling laws deduced by iDe Beck et alj 
I2OI0I) . 





Figure 7. HCN observed spectral lines (gray) are compared to 
the spectral line predictions (black) based on the CSE model 
shown in Fig. |4] and the abundance stratification displayed in 
Fig.m 



Comparison to other observational studies: Most observa- 
tional studies listed in Table|6]assume a fractional CO abundance 
of [CO/H] = 1 - 1.5 X 10"''. Using different CO rotational lines 
the mass-loss rate is then derived. The deduced ^^CO/^'^CO ratio 
ranges between 9 and 14. 

Using CO rotational line transitions, other studies have also 
estimated the mass-loss rate (see Table |7]i. The results depend 
on the assumed or derived temperature distribution, the dis- 
tance, the adopted [CO/H2] abundance ratio, and the radiative 
transfer model or analytical approximation used. All (scaled) 
mass-loss rate values are in the narrow range between 6.5 x 
10~^ and 9 x 10~^ Mp/yr, th e exception being the result of 
iGonzalez Delgado et al.l (l2003h . which is a factor ^4 higher We 
note that the work of these authors was not devoted to the study 
of CO, and it remains unclear what line intensities they used in 
their modeling. 

4.3.2. HCN 

Results: The narrow Gaussian line profiles clearly point to- 
ward (at least) an inner wind origin of HCN. As discussed above, 
we impose a value R2 = 500 i?* to simulate the mapping re- 
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Table 6. Comparison of the deduced fractional abundances to other observational studies and theoretical predictions. All fractional 
abundances are given relative to the total H-content. In cases where values found in literature were given relative to H2, they were 
re-scaled relative to the total H-content by assuming that all hydrogen is in its molecular form H2. 







*^C0 


cs 


HCN 


CN 


SiO 


SiS 


SO 


SO2 


Lindavist et al. (1988)" 
Omont et al. (1993)" 
Buianabal" al. (1994)"= 
Kim et al. (2010)" Cease A") 
Kimetal. (2010)"^ Cease B') 


1.5 X 10"* 

1.5 X 10"* 
1.5 X 10""* 
1.5 X 10~* 


1.6 X 10"^ 
1.75 X 10"^ 
1.75 X 10"^ 


1.5 X 10"'' 

5.0 X 10"** 

3.0 X 10"'' 

8.1 X 10"** 


3.0 X 10"'' 

4.9 X 10"^ 
1.4 X 10"'^ 
4.3 X 10"'' 


- 

1.6 X 10"^ 
5.1 X 10"** 


1.5 X 10"" 
8.5 X 10"" 
1.3 X 10"^ 
5.1 X 10"" 


3.5 X 10"'' 

2.2 X 10"^ 

1.3 X 10"" 
3.1 X 10"^ 


9 X 10"' 
1.3 X 10"" 
7.8 X 10"'' 
2.7 X 10"^ 


2.05 X 10"" 

1.4 X 10"^ 
4.2 X 10"" 


Gonzalez Delgado*^ 

et al. (2003) 

Schoiei- et al. (2007a) ^ 

this work 


1.0 X 10"* 
1.0 X 10"* 

1.0 X 10"* 


7.1 X lO"'' 


4 X 10"** 


2.2 X 10"'' 


1.0 X 10"^" 
-3.0 X 10"** 


2.0 X 10"" 

8.0 X 10"" 
-2.0 X 10"'' 


5 X 10"" 
-5.0 X 10"" 
5.5 X 10"" 
-4.0 X 10"** 


2.0 X 10"^ 


1.0 X 10"" 


Duariet al. (1999)'' 
Chei-chneff (2006)' 
Wmacv& Millar (1997)^ 


5.38 X 10"* 
6.71 X 10"* 
4 X 10"* 




2.75 X 10"^ 
1.85 X 10"^ 
2.9 X 10"''' 


2.12 X lO"'' 
9.06 X 10"^ 
1.4 X 10"'' 


2.40 X 10"^° 
3 X 10"**' 
3.5 X 10"'' 


3.75 X 10"^ 
4.80 X 10"^ 
3.2 X 10"^ 


3.82 X 10"*° 

7 X 10"** 
3.5 X 10"" 


7.79 X 10"** 

1 X 10"^ 
9.1 X 10"'' 


1 X 10"*^ 
2.2 X 10"" 



Notes. In the first part of the table, observational resuUs are listed based on the assumption of optically thin emission and a population distribution 
which is thermalized at one excitation temperature. The second part gives observational r esults based on a n on-LTE radiative transfe r analysis. 
Theoretical predictions for either the inner envelope jPuari et"ai 1 119991 : ICherchneffll200^ or outer envelope j'Willacv & Millajll997h fractional 
abundances are given in the last part. 

References. 

no information on used distance and mass-loss rate 
distance is 270pc, M = 4.5 x 10"® Mq/yi 
distance is 270 pc, Af = 4.5 x 10"*^ M^/yi 

distance is 250 pc, assumed Al of 4.7 x 10"® Mq/jt, LTE is assumed, 'case B' represents a solution with a larger outer radius than 'case A' 

Te in Gaussian distribution for SiO is 2.5 x 10^® cm, distance is 250 pc and M = 3 x 10"'' Mq/ji 
^'.rein Gaussian distribution for SiS, distance is 260 pc and Af = 1 x 10"^ Mq/yi. For 2-component model: fcis5.5 x 10"" and taken constant 
out to 1.0 X 10^^ cm and the lower abundance Gaussian component has fo of 5.0 x 10"^ and re of 1.6 x 10^® cm. Using one (Gaussian) 
component distribution, fo is 5 x 10"* and re is 1.6 x 10^® cm. 
^ : only value at 5 7?* is given 

'': predicted values at 2.2 i?* in the envelope for IK Tau 
': predicted values at 2 in the envelope for TX Cam 
: predicted peak fractional abundances in the outer envelope 



suits of iMarvell (l2005h . The deduced abundance around 50 i?* 
is [HCN/H] = 2 - 2.5 X 10"^. 



Compari son to theoretica l pre dictions: Theoretic al predic- 
tions by iDuari et al ] (Il999h and [ Cherchneffl (l2006l) and ob- 
servational studies by, e.g., iBieging et al.l ( 200ol) and iMarvell 
(l200l indicate that HCN forms in the inner wind region of 
M-type envelopes by shock-induced non-equilibrium chemical 
processes. This is cont rary to the photochemical models of 
IWillacv & Millaii (|1997|) . where HCN is only produced in the 
outer envelope by photochemical reactions involving NH3 and 
CH3, which they assumed to be parent species originating close 
to the stellar photosphere and injected in the outer region. This 
formation route, however, leads to line shapes which are clearly 
non-Gaussian but are flat-topped (for optically thin unresolved 
emission) since the wind is already at its full velocity in this re- 
gion. 

The maximum size distribution of 3.85" derived bylM arvell 
(l200l is suggested to be caused by photodissociation of HCN. 
iBieging et al. (l2000h used a parametrized formula to describe 
the photodissociation radius for HCN as a function of gas mass- 
loss rate and wind velocity (see their Eq. 2). This estimate leads 
to a HCN photodissociation radius of 1.2 x 10^^ cm, in good 



agreement with the result o f|M arve il (l2005h of 7.6 X 10 cm (at 
265 pc). Using HCN as a parent species with an injected abun- 
dance of 1.5 X 10"^, new chemical outer wind models were cal- 
culated using the code as described in , Willacy & Millar (199^. 
The derived photodissocation radius is around 2 x 10^^ cm, be ing 
a factor ^3 higher than the observed value o f|M arve l (l200l . 

Compared to the t heoretical predictions for TX Cam at 5 
by ICher(:hnef j (l2006h or for IK Tau at 2.2 by iDuari et al.l 
(I1999I) . our deduced abundance is a factor 10 to 40 lower, re- 
spectively. There are a few possibilities for th e origin of this dif - 
ference. (i.) Contrary to what is thought (e.g. lDuari et al.|[T999l) . 
HCN may participate in the formation of dust grains in the inner 
envelope, (ii.) The formation mechanism of HCN in the inner 
wind is directly linked to its radical CN by 



CN + H2 ^HCN + H. 



(4) 



The destruction route is the reverse reacti on. Shocks trigge r 
CN and further HCN formation in the gas (ICherchneffll2006h . 
The formation processes of both molecules depend critically on 
the physical parameters of the shocked gas, specifically on the 
physics in the 'very fast chemistry zone', that is the naiTow re- 
gion after the shock front itself. The modeling of this zone is 
still subject to many uncertainties (e.g. cooling rate, velocity. 
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12C160(1-0) (OSO) 12C16O(1-0) (IRAM) 12C160(2-') (JCMT) 12C160(2-1) (IRAM) 




20 40 60 80 20 40 60 80 20 40 60 80 20 40 60 80 

Velocity [km/s] Velocity [km/s] Velocity [km/s] Velocity [km/s] 



12C160(3-2) (JCMT) 1 2C1 60(3-2) (APEX) 1 2C' 60(3-2) (APEX) 1 20 1 60(4-3) (JOMT) 




20 40 60 80 20 40 60 80 20 40 60 80 20 40 60 80 

Velocity [km/s] Velocity [km/s] Velocity [km/s] Velocity [km/s] 



120160(4-3) (APEX) 12016O(/-6) (APEX) 130160(2-^) (JOMT) 130160(3-2) (APEX) 
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Figure 6. CO rotational line profiles of IK Tau (plotted in grey) are compared with the GASTRoNOoIVI non-LTE line predictions 
(in black) using the parameters of 'model 3' as specified in Table |4] The rest frame of the velocity scale is the local standard of rest 
(LSR) velocity. 



shock strength), yielding uncertainties on the theoretical frac- 
tional abundances of at least one order of a magnitude. 

Comparison to observational studies: All observational de- 
duced values agree within a factor ^2, but are clearly lower than 
the inner wind non-TE theoretical predictions (see Table|6l). 

4.3.3. SiS 

Results: Next to CO, SiS is the molecule with the most molec- 
ular line transitions at our disposal. From the SiS(5^) up to 
SiS(20-19), excitation temperatures from 13 to 183 K are cov- 
ered, and one can trace the envelope abundance spatial variations 
between 40 and 6000 i?*. The lower lying transitions (SiS(5^) 
and (6-4)) probe material at larger radial distances from the star; 
the SiS(19-18) and (20-19) line intensities are sensitive to the 
abundance in the inner wind region. IVIodeling the strength of 
the low and high-excitation lines gives evidence of a depletion of 
SiS: the abundance at 40 R-^- is estimated to be 5.5 x 10^® and 
decreases to ^ 4 x 10^^ at 120 i?*. This is also evidenced by the 
extra sub-structure in the high signal-to-noise high-resolution 



line profiles of the higher excitation lines observed with JCMT 
and IRAJVI, where the small peak traces a high abundance com- 
ponent at low velocity close to the star and the broader plateau 
the lower abundance component further away in the envelope. 
SiS can be depleted due to the adsorption of SiS molecules onto 
dust grains, before being photodissociated much farther away. 
Using ISI 1 1 yum interferometric data. Hale et al. ( 1997) detected 
dust excess, with a first intensity peak around 0.1", going out 
up to at least 0.7"(or 2.8 x IQi^cmw 180 i?*) around IK Tau. 
SiS condensation onto dust sp ecies has been postulated already 
(e.g. iBieging&Tafallall 19931) . but not all atomic Si or molecu- 
larSiS should condense onto dust species (e.g. iDominik et al.l 
Il993h . The different rotational lines give us the potential to con- 
strain the compact, pre-condensation, SiS fractional abundance 
component quite well, but the uncertainty on the extended post- 
condensation component is higher. 



Comparison to theoretical predictions: The abundance at 
40 i?* is higher than the theoret i cal va lue predicted in the in - 
ner wind envelope bv lDuari et"an d 19991) and lCherchnefj ( l2006l) . 
This may either indicate that the destruction reaction which oc- 
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Table 7. Mass-loss rate values derived from ^^CO rotational line transitions for IK Tau. The first and second columns list the 
distance and [CO/H2] value used in the different studies, the third column gives the derived gas mass-loss rate, the fourth column 
contains the mass-loss rate scaled to our adopted value s for the distance (D = 265pc) and CO-abundance ([CO/H2] = 2 x lO^'*) 
using the scaling laws deduced bv lDe Beck et al.l (1201 Ol) . the fifth column lists the rotational Une transitions used, and last column 
gives the reference. 



D 


[CO/H2] 


M 




M (scaled) 


Lines 


Reference 


[pc] 






[Mo/yr] 




[Mo/yr] 






350 


2 X 10" 


4 


1.0 X 10" 


5 


6.7 X lO"*^ 


^^CO(l-O) 


Buiarrabal et al. (1989) 


270 


5 X 10" 


4 


4.0 X 10" 


6 


6.8 X lO"*^ 


^^CO(l-O) 


SoDkaetal. fl989) 


260 


5 X 10" 


4 


4.4 X 10" 


6 


7.9 X 10"^ 


^^CO(l-O) 


Loucetal. (1993) 


260 


5 X 10" 


4 


3.8 X 10" 


6 


6.8 X 10"^ 


^^CO(I-O) 


Neri et al. (1998) 


250 


2 X 10" 


4 


3.0 X 10" 


5 


3.2 X 10"^ 


^2 CO( 1-0)^(4-3) 


Gonzalez Delaado et al. (2003) 


250 


3 X 10" 


4 


4.7 X 10" 


6 


6.5 X lO"*^ 


^2cO(l-0)->(3-2) 


Tevssier et al. (2006) 


300 


2 X 10" 


4 


1.0 X 10" 


5 


7.7 X lO"*^ 


^2cO(l-0)->(4-3) 


Ramstedt et al. (2008) 


265 


2 X 10" 


4 


9.0 X 10" 


6 


9.0 X 10"^ 


^2CO(l-0)->(7-6) 


this study 



curs a t larger radii in the inner O-rich envelope via dCherchnefj 
l2006h 



S + SiS ^ S2 + Si 



(5) 



is not as efficient, or that an extra formation route is not yet taken 
into account in the theoretical modeling. In general, the uncer- 
tainties on chemical reaction rates involving sulfur are still very 
high (I. Cherchneff, priv. comm.). The high abundance around 
AORj, is in accordance with the observational suggestion by 
iDecin et"al ] f2008b) that SiS forms close to the star, whatever 
the C/O ratio of the target. 

Comparison to observational studies: Similar results con- 
cerning the depletion of S i S in th e intermediate wind region are 
obtained bv lSchoier et al.l ( l2007al) . From the simplified analyses 
assuming optically thin emission thermalized at one excitation 
temperature (see Table |6]l it is not possible to derive this kind of 
abundance depletion pattern. 

4.3.4. SiO 
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Figure 10. Upper pan e l: ^^S iO spectra of IK Tau from 
iGonzalez Delgado et al.l (l2003l) . Lower panel: ^^SiO spectral 
line predictions based on the CSE model shown in Fig. |4] and 
the abundance stratification displayed in Fig.|5] 



^^SiO(7-6) - APEX 
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Figure 11. ^^SiO observed spectral lines (gray) are compared to 

naso-^i. J .11- / \ J. the spectral line predictions based on the CSE model shown in 

i igure 9. °SiO observed spectral lines (gray) are compared to r- m 1 1 1 1 ■ i- 1 1 • i-- m 

.II- J- u J /-cV J 1 u Fig. MJ and the abundance stratification displayed in Fig. dI 

the spectral line predictions based on the CSE model shown in t- j a <—i 

Fig.|4]and the abundance stratification displayed in Fig.|5] 



Results: For the ^^SiO isotopolog, 5 transitions were observed, 
with excitation energies ranging between 6 and 75 K. The line 



strengths and profile shapes of all ^^SiO, ^^SiO, and •^°SiO lines 
are well predicte d, exce pt for t he ^^SiO(6-5) lin e as observed 
with the SEST bv lGonza lez De lgado et al.l (l2003h (see Fig. [TOll. 
Since the strength of both the 2^SiO(5^) and 28SiO(7-6) are 
well reproduced, and both of these lines share the line formation 
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Figures. SiS observed spectral lines (gray) are compared to the spectral line predictions based on the CSE model shown in Fig.|4] 
and the abundance stratification displayed in Fig. [5] 
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Figure 12. '^''SiO observed spectral lines (gray) are compared to 
the spectral line predictions based on the CSE model shown in 
Fig. |4] and the abundance stratification displayed in Fig.|5] 



region with the ^*SiO(6-5) transition, an absolute calibration un- 
certainty can be the cause of this discrepancy, but time variability 
of the lines should also be considered (see Sect. 15. 11 1. 

Although not as pronounced as for SiS, the modeling of the 
different rotational transitions indicates an abundance decrease 



with a factor ^40 around 180 i?*. SiO is a parent molecule and a 
strong candidate to be depleted in the wind of O-rich envelopes: 
at la rger radii in the in ner envelope, OH alters the SiO abundance 
via (ICherchneffll2006h 



SiO + OH ^ SiOa + H . 



(6) 



Si02 may condense as silica. It may also participate in the for- 
mation of amorphous or crystalline silicates. 

The derived isotopic ratios in the envelope are 
[28SiO/2^SiO]~27 and [^SSiO/^^SiO] = 80, [^^SiO/^OSiO] = 3. 
They are discussed in Sect. 15.21 

Comparison to theoretical predictions: The SiO abundance 
at the inner dust condensation radius is slightly below the theo- 
retical predictions of Duari et al. (1999) and Cherchneff (200^, 
which is very reasonable taking the assumptions of both the the- 
oretical chemical kinetic calculations and our modeling into ac- 
count. It possibly points toward the condensation of SiO onto 
dust grains in the intermediate wind zone, before 70 i?*, a re- 
gion where our observed lines are not very sensitive to the exact 
abundance distribution. The observational study bv lDecin et alj 
(I2008bh indicates that SiO is formed close to the star, in sup- 
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port of the theore tical predictions by iDuari et al ] (Il999h and 
ICherchnefB (l2006l) . In the outer envelope, SiO is very stable 
and only photodissociat ed around a few thousand stellar radii 
(IWillacv&Millailll997l) . 



Comparison to observational studies: iLucas et all d 19921) 

mapped the ^^SiO(2-l) u = flux distribution, showing that it 
has a circular geometry. The half-peak intensity radius has a di- 
ameter of 2.2 +1- 0.1" (or a radius of 4.35 x 10^^ cm at 265 pc, 
in our model being 190 i?*). The SiO(2-l) emission regions in- 
dicate that the final expansion velocity is not yet reached in the 
SiO emission region, suggesting that grain formation must still 
take place as far as 10^® cm from the star This supports our re- 
sults on the velocity structure based on the study of the HCN line 
profiles (Sect.l4Tli. 

Compared to other observational studies, the deduced abun- 
dance around 70 Ri, is quite h igh, while the outer wind abun - 
dance agrees with the result bv lGonzalez Delgado et al.l (l2003l) . 



to model (see Sect. 14.3.2b and that the sulfur reaction rates are 
not well known (see Sect. 14.3.3b . this difference is not so cum- 
bersome. However, as suggested for HCN, it may also be the 
case that CS or the radical CN are involved in dust formation, 
altering its abundance in the intermediate wind region. The low- 
resolution of the two observed CS lines do not provide enough 
information to firmly prove this. In the outer envelope, CS is 
first formed from H2S. Somewhat farther away, the reaction of 
atomic carbon with SO a nd HS forms CS, before it is photodis- 
sociated by UV radiation dWillacv & Millaijll997h . 



Comparison to observational studies: While our deduced 
CS fractional abundance agrees with the values derived by 
.Bujarrabal et al. ( 19 94) and K im et al., C201(3i) . the observational 
result bv lLindqvist er'al.| (1 1 988 ) is higher by a factor ~4. 



4.3.6. CN 



4.3.5. CS 
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Figure 14. CN observed spectral lines (gray) are compared to 
the spectral line predictions based on the CSE model shown in 
Fig. m and the abundance stratification displayed in Fig.|5] The 
Figure 13. CS observed spectral Unes (gray) are compared to flie dashed line predictions corresponds to the 'alternative solution' 
spectral line predictions based on the CSE model shown in Fig.|4] shown in Fig. [5] 
and the abundance stratification displayed in Fig.|5] 



Results: For CS, we only have two lines at our disposal, the (6- 
5) and (7-6) rotational transitions, both with a low S/N-ratio. The 
fractional abundance at 300 i?* is estimated to be 4 x 10^^. 

Comparison to theoretical predictions: The derived abun- 
dance of ^ 4 X 10^^ is much higher than the TE-abund ance 
of 2.5 X 10~^^ for an oxygen-rich star ICherchnef j (l2006l) pre- 
dicts CS to be a parent molecule, with a non-TE abundance for 
TX Cam around 1.8 x IQ-^ at I R^ and around 2 x lO"'^ at 
5 i?* away in the envelope, while iDuari et aP d 19991) predicts a 
value of 2.75 x 10^^ at 2.2 Ri,. In thei r theoretical modeli ng of 
the carbon-rich AGB-star IRC-Hl0216 lMillar et alj (1200 lb also 
argue for the need of CS as a parent species to account for the 
vibrationally excited CS lines detected in IRCh-10216. 

As in the case of HCN (Sect. 14.3.2b . our derived abun- 
da nces are a f a ctor ^ 50 lower compared to the predictions 
of 'Cherchneff (200^ and a factor ^8 lower compared to 
[Pua ri et al. (1999). The dominant formation pathways of both 
CS and HCN occur in the fast chemistry zone of the gas parcel 
excursion involving CN. Knowing that this zone is very difficult 



Results: The CN lines display a peculiar profile, probably re- 
lated to the hyperfine structure of the molecule. Alt hough the 
signal -to-noise of the individual components is low, iKim et aP 
(120101) noted already that the strength of the different peaks are 
not in agreement with the optical thin ratio of the different hy- 
perfine structure co mpoments and hint to h yperfine anomalies as 
already reported bv lBachiller et al.l (Il997h . Simulations with the 
GASTRoNoOM-code taking all the hyperfine components into 
account confirm this result. We therefore opt to simulate both 
CN lines using the strongest component only. I.e., for the N = 
3-2, J ^ 5/2- 3/2 line we will use the F ^7/2- 5/2 com- 
ponent at 34003 1.5494 MHz, for the AT = 3-2, J = 7/2-5/2 
linetheF = 9/2-7/2component at 340248.5440 MHz is used. 

Due to low signal-to-noise ratio of the lines and the prob- 
lems with the different hyperfine componets, the derived abun- 
dance fractions are loosely constrained. To illustrate this, two 
model predictions are shown in Fig. [141 For one model, the in- 
ner abundance ratio is taken to be 3 x 10~® and from 1000 i?* 
onward, the abundance stratification follows the predictions by 
Willacy &_Millar (1997) (dashed line in Fig. |5]l. For the other 
model, the inner abundance ratio, /i, is lowered to 1 x 10^^" 
yielding a peak fractional abundance around 2000 i?,^ of 2.5 x 
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10 ^, being a fac tor ~8 higher tha n the p eak fractional abun- 
dance derived by IWillacv & Millar! (Il997h (see dotted line in 
Fig.|5). 



crease the predic t ed pe ak fractional abundance computed by 
IWillacv & Miii^ (Il997h which was somewhat too low com- 
pared to the observed value hsted in their Table 6. 



Comparison to theoretical predictions: The predicted inner 
wind abundance fractions give higher pref erence to the sec- 
ond model described in previous paragraph: iDuari et al ] (ll999h 
predicts a fractional abundance of CN around 2.4 x 10"^" at 
2.2 Ri, for IK Tau, while Cherchneff (2006) gives a value of 
3 X 10-1" at 5 i?* for TX Cam. Further out in the wind, CN is 
produced by neutral-neutral reactions and the photodissociation 
of HCN, yielding peak fractional abundances around 3 x 10^^ 
(IWillacv & Millaiill997h . 

If HCN is indeed photodissociated around 400 to 500 i?* 
(see Sect. I4.3.2l i. the peak fractional abundance of CN is not ex- 
pected to occur around 2000 i?* (see the model predictions by 
IWmacv&Miiia3ll997l in Fig.©, but shghdy beyond 500 i?,^ 
since the photodissociation of HCN is the main formation route 
to CN in the outer envelope. Shifting the CN peak fractional 
abundance in the second model to 500 i?* with an abundance 
value of 1.5 x 10"'^ also yields a good fit to the (noisy) data. 

Comparison to observational studies: This is the first time 
that the non-LTE CN fractional abundance for IK Tau has been 
derived, although the uncertainty on the derived abundance is 
large. 



Comparison to observational studies: For the first time, the 
SO abundance fraction is derived using a non-LTE radiative 
transfer analysis, although the low S/N prevents an accurate 
abundance determination. The result agrees with the LTE anal- 
ysis by Kini et al. ( 20i3L but is a factor of a few lower than 
lOmont et al.l (tl993b and lBuiarrabal et alj (|1994|) . 

4.3.8. SO2 

As for SO2 the main formation channel in both the inner and 
outer wind region is 

SO + OH ^ SO2 + H . (7) 

In the outer wind , SO? rapidly photodissociates back to SO 
(IWillacv & Millad |1997[) . The inner wind predictions for TX 
Cam only yield an abundance of 1.3 x 10" but the detec- 
tion o f SO2 at 7.3 /im in several O-rich giants dYamamura et al.l 
1 19991) would imply a formation site close to the star and an abun- 
dance in the range IQ-* — 10"^, slightly below the SO values. 
ICherchneffI (l2006l) argues that a limited number of reactions in 
the involved SO2 formation scheme and the lack of measured 
reaction rates may explain this discrepancy. 



4.3.7. SO 
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Figure 15. SO observed spectral lines (gray) are compared to the 
spectral line predictions based on the CSE model shown in Fig.|4] 
and the abundance stratification displayed in Fig.jS] 



Results: As in case of CN (Sect. 14.3.6b the noisy profiles pre- 
vent a good determination of the abundance stratification. A 
fractional abundance of ^ 2 x 10^^ is derived around 200 i?*. 

Comparison to theoretical predictions: An abundance 
stratificat ion compat i ble wi th both the inner wind predic- 
tions of ICherchnefj (l2006h and the outer wind model of 
IWillacv & MiiiaT " (1 19971) can be derived yielding a good 

representation of both SO lines observed with APEX. 

fWillacv&MiTia3(ll997h assumed no SO injection from the in- 
ner wind at large radii, but in-situ formation processes only. 
Assuming that SO is indeed injected to larger radii can in- 



Results: The availability of 10 different transitions gives some 
hope that we can shed light on the discussion about the in- 
ner wind fractional abundance. However, it turns out that we 
are unable to fit the high-excitation S02(17i 17 — 16o,i6) and 
S02(132.i2 — 12i,ii) lines observed with APEX (see below). The 
high-excitation S02(144^io - 143^0 and S02(1434i - 142, 12) 
can be predicted quite well. These two lines are narrower than 
the other SO2 lines in the sample, indicating that their full for- 
mation region is in the inner wind region where the wind has not 
yet reached its full expansion velocity. 

Extensive modeling efforts were done to predict the high- 
excitation SO2 observed with APEX. While the high-excitation 
lines involving the J = 14 — 14 levels are reasonably well pre- 
dicted with the models proposed above, the J =13 — 12 and 
J = 17 — 16 are far too weak. In a study of SO2 in star form- 
ing regions, fvan der Tak et al ] (120031) encountered an analogous 
problem, which they solved by introducing a high temperature, 
high-abundance component. The increase in abundance could be 
a factor of 100-1000. Introducing an unrealistic compact, very 
high-abundance component with /i > 1 x 10^^ up to 200 Ri, 
reproduces the APEX J = 13 — 12 and J = 17—16 within 
a factor 2, but the J = 14 — 14 lines, involving similar ex- 
citation levels, is a factor -^15 too strong. A possible cause of 
the discrepancy could be a misidentification of the observed 
lines. However, different line data bases always point towards 
an identification as SO2 transitions. But also the collision rates 
from the LAMDA database ma y be problematic. In the LAMDA 
database, ISchoier et aL ("2005") have extrapolated the collisional 
rates as computed by|Green (1995). Green ( 1995^ computed the 
collisional rates for temperatures in th e range from 25 to 12 5 K 
including energy levels up to 62 cm~ ^ : ISchoier et aH (l2005l) ex- 
trapolated this set of collisional rates to include energy levels up 
to 250 cm^i and for a range of temperatures from 10 to 375 K. 

Neglecting the S02(17i,i7-16o,i6) and S02(132,i2-12i,ii) 
lines, the strength of the other eight SO2 lines can only be ex- 
plained using a high inner abundance ratio of 1 x 10"^, clearly 
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Figure 16. SO2 observed spectral lines (gray) are compared to the spectral line predictions based on the CSE model shown in Fig. |4] 
and the abundance stratification displayed in Fig. [5] 



poin ting toward an inner w ind formation region (in accordance 
with lYamamura et al.ll 19991) . 



Comparison to observational studies: As for SO, the SO2 
fractional abundance is derived for the first time using a full non- 
LTE radiative transfer analysis. The derived abundance value is 
somewhat lower than the results from classical studies assuming 
optically thin emission and one excitation temperature, although 
we have to state clearly that the SO2 modeling still poses many 
problems. 



5. Discussion 

The derived fractional abundances are already discussed in 
Sect. 14.31 In this section, we focus on the possible time variabil- 
ity of the emission lines and on the derived SiO isotopic ratios. In 
the last part, H2O line profile predictions for the Herschel/HIFI 
mission are performed. 



5.1. Time variability 

The observed molecular emission lines could be time variable. 
Unfortunately, no dedicated study has yet been performed to 
study the time variab ility of the molecular lines in IK Tau. 
ICarlstrom et alj d 19901) reported on a monitoring programme of 
the SiS(J = 4-3, 5^, and 6-5) emission from the Mira-type 
carbon-rich AGB star IRC H- 102 16. It was found that the circum- 
stellar J = 5-4 and J = 6-5 line emission toward IRC H- 102 16 
varies both in line intensity and in line shape. A clear correlation 
between the variations and the infrared flux (as measured us- 
ing the /■C-band magnitude) is found for the J = 5-4 and J = 6-5 
lines (Bieging & Tafalla 1993), but not for the J = 4-3 line. This 
indicates that, at least, the population of a few levels vary in 
phase with the stellar flux. A change in the pumping mechanism 
of the infrared vibrationally excited levels of a molecule will 
modify the excitation in the ground vibrational state. A change in 
dust emission may also alter the excitation of a molecule, since 
dust emission has the potential of affecting the level populations 
in the ground vibrational state, in particular through pumping via 
excited vibrational states. 

ICernicharo et al.l ( l2000h have looked for time-related inten- 
sity variations in a line survey at 2 mm of IRCh-10216. Among 
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the 2-mm hnes, the most likely lines to be affected are: (i) CS, 
HC3N, SiO and SiS, four species whose IR lines are known to 
be optically thick, as well as (ii) the vibrationally excited lines 
of C4H and HCN. During the 10-year long run, these lines were 
observed at several occasions. The ground-state mm lines were 
found to have stable shapes and intensities (within 20 % which 
is consistent with the calibration uncertainty). Only the strong 
HCN, 1/2 — I, J = 2 — 1 line, which is known to be masering, 
showed a factor of 2 intensity variation with time. 

Currently the effect of time variability on circumstellar line 
emission is unknown for AGB stars in general. 

5.2. SiO isotopic ratios 

The SiO abundance isotopic ratios derived for IK Tau are 
28SiO/29SiO = 27, 28siOPSiO = 80, and 29SiO/30SiO = 3, with 
an uncertainty of a factor of ^2 due to the low signal-to-noise 
ratio of the ^SSiO and ^^SiO lines. The 29SiO/''^°SiO is similar 
to the simple ratio of integrated intensities corrected for a dif- 
ference in transition strength and beam filling factor (the com- 
bined effect is a frequency factor of i/^^), but the ^^SiO/^^SiO 
and ^^SiO/'^'^SiO ratios are a factor ^1 and ^10 larger, respec- 
tively, due to neglect of optical depth effects of the ^^SiO lines in 
case the simplified intensity ratio is used. Compared to the solar 
isotopic ratios of (^^SiO/^^SiO)© = 19.6, (28SiO/3OSiO)0 =29.8 
and (29SiO/^°SiO)o = 1.52, IK Tau is underabundant in neutron- 
rich isotopologs or overabundant in •^^Si. 

Due to the weakness of the ^^SiO and ^°SiO lines, only 
few results on the silicon isotopic ratios in the circumstellar 
envelopes around AGB stars are reported in lite rature. By fit- 
ting the SiO maser emissivity Ijiang et aP (l2003h estimated the 
^^SiO/'^^SiO isotopic ratio in two oxygen-rich Miras, R Cas 
£ind N V Aur, as being 29 and 32, respectively. lUkita & Kaiful 
(Il988h measured the relative intensities of the ^^SiO/'^"SiO lines 
(J = 2 - l,v = 0), being 2.4 for the S-type Mira x Cyg 
(with C/O ratio slightly lower than 1), 1.5 for IK Tau, and 2.9 
for the oxygen-ric h Mira Villi Op h. For the carbon-rich AGB 
star IRCh- 10216. ICernicharo et all ( l2000h and iHe et alj ( l2008h 
derived respectively 28SiO/29SiO= 15.4 ± 1.1 (17.2 ± 1.1) and 
^^SiOPSiO = 20.3 ± 2.0 (24.7± 1.8) and ^SSiO/^OSiO = 1.45 ± 
0.13 (1.46 ± 0.11) from the ratios of integrated line intensities, 
being close to the solar values. The results on IRCh-10216 are, 
however, lower limits, since no correction for opacity effects has 
been done. 

iLambert et alj ( Il987h analysed high-resolution spectra of the 
SiO first overtone band around 4 /im. They obtained estimates 
of the atmospheric ^^SiO/^^SiO abundance ratios for four red 
giants. For the M-type /3 Peg and the S-type star HR 1105, the 
^^SiO/^^SiO ratio is close to the solar ratio. ^^SiO appears to 
be underabundant in the MS star Ori (2^SiO/2^SiO = 40) and 
the M-type star 10 Dra (2*^SiO/29SiO~53). The ^^SiO isotope 
appears to be underabundant by a factor of ^2 in all four red 
giants. 

iTsuii et aP (1 19941) reported on high spectral resolution ob- 
servations of the 4 /im SiO first overtone band in six late-type 
M giants and two M supergiants. The atmospheric •^^Si/^^Si 
and ^^Si/'^'^Si and ^^Si/^^^Si ratios in the M giants are always 
slightly lower than the terrestrial values, i.e. more neutron-rich 
nuclei tend to be more abu ndant. This is opposite to the re- 
sults of iLambert et al.l (Il987h . and assuming that the isotopic ra- 
tios are not modifi ed in the circumstellar envelope, the result of 
iTsuii et aTl d 19941) is also in contrast to the circumstellar isotopic 
ratios of M-type giants listed above. 



The silicon isotopic ratios are not obviously correlated with 
other stellar properties. In the literature, it is conventional to ex- 
press the silicon (and other element) isotope ratios in parts per 
thousand deviation from the solar silicon isotope ratio: 
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yielding values for IK Tau of —274 and —627, respectively. 
These values reflect both the initial isotopic composition of the 
star and possible effects due to nucleosynthesis. 

For an AGB star, the initial Si isotopic composition in the en- 
velope is altered by slow neutron capture reactions (s-process) 
in the He intershell and subsequent third dredge-up (TDU) 
events, increasing the ^^Si and ■^"Si abundance fractions. A low- 
metallicity star is expected to have initial ^^Si/^^Si and ''"Si/^'^Si 
ratios that are smaller than the solar ratios. The inferred isotopic 
shifts of the Si isotopes are smaller for an O-rich than for a C- 
rich AGB star since a C-rich star goes through mo re dredge-up 
events increasing th e ^^C and s-processed material dZinner et al.l 
2006; VoUme r et aT] 2008 ). Using two different stellar evolution 



codes. IZinner et al.l (l2006h studied the change in silicon isotopic 



ratios in AGB stars: the shift in Si isotopic ratios and the increase 
of the ^^C/^^C in the envelope during third dredge-up are higher 
for higher stellar mass, lower metallicity, and lower mass-loss 
rate, but their predicted silicon isotopic shifts are always much 
higher than the observational values derived for IK Tau. The 
mimimum values plotted in their Fig. 6 correspond to the value 
(5© Si) and (5o(^°Si) at C/0=1, and are larger than -200. They 
find that no noticeable changes in the Si isotopes occur when the 
star is still O-rich. Consequently, the isotopic anomalies in sili- 
con found for several M-type giants probably reflect those of the 
interstellar medium out of which stars were formed. 

For the solar system material, the silicon isotopic ra- 
tios are thought to be understood by a mixture of the nu- 
clear products by types I and II supernovae dTsuii et al.l 
II994I) . The predicted silicon isotopic ratios by type II super- 
novae are around 2*Si/29Si~15 an d 2^Si/^°S i ^35 according 
to JHoppe et al. (2009), while Hashi moto et al.l (Il989h arrive at 
lower values being 8.9 and 12.6 respectively. Type I supernovae 
produce mostly ^^Si with Uttle ^^Si and ■^"Si ( [Thielemann et al.l 
|T98a). Non-terrestrial silicon isotopic ratios can then be reason- 
ably explained in the same way as for the solar system but by 
assuming a different contribution of types I and II supernovae. 

The silicon isotopic shifts reported here for IK Tau are 
much lower than v alues de duced from presolar silicate grains 
(I VoUmer et al.ll200 8; Moste faoui & Ho ppe 2004), of which the 
origin spans the range from red giant branch (RGB) and AGB 
stars up to supernovae. Looking to s ilicon isotopic ratio s derived 
from presolar SiC grains (Fig. 2 in IZinner et al.ll2006l) . the sili- 
con isotope ratios of IK Tau correspond to the X-grains, which 
are thought to originate in Type II supernovae (however major 
discrepancies be tween model predi ctions and observed isotopic 
ratios still exist; iNittler et al1ll995h . The ^^C/^^C ratio inferred 
for IK Tau (=14) is at the l ower Umit of the v alues derived for 
X-type grains (see Fig. 1 in lZinner et al.ll20d6h . 

Hence, if the atmospheric (and circumstellar) silicon isotopic 
ratio is indeed not changed due to nucleosynthesis and subse- 
quent dredge-ups, the above arguments seem to suggest that the 
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interstellar medium out of which IK Tau was born has a mix- 
ture analogous to X-type grains of which supernovae type II are 
thought to be the main contributors. The measurement of other 
isotopic ratios can shed new light on this discussion. 



5.3. H2O line profile predictions 

Line profile predictions are performed for a few water lines 
which will be observed by Herschel/HIFI in the framework 
of the Guaranteed Time Key Programme HIFISTARS (P.I. V. 
Bujarrabal) (see Table |8]i. This key programme focusses on the 
observations of CO, H2O and HCN lines in a well-selected 
sample of evolved stars in order to gain deeper insight into 
the structure, thermodynamics, kinematics and chemistry of 
CSEs and into the mass-loss history of evolved stars. The in- 
ner wi nd abun dance fraction is assumed to be [H20/H2]=3.5 x 
10~* (ICherchnefti2006). The photodissociation radius is taken 
from the modeling of Willacv & MillaH (Il997h . being around 
1600i?^. Using the analytical fo rmula from &oenewe genl 
( 119941) deduced from the results of iNetzer & Kna ppI (j 19871), a 
photodissociation radius of 2.8 x 10^^ cm or 1870 i?* would 
be obtained. As standard set-up, the Barber H2O line list (see 
AppendixlAb is used, including 45 levels in both the ground state 
and first excited vibrational state (the bending mode 1^2 = 1 at 
6.3 ^m). 



H2O{1,„-1„,),v = 0-0 



H,0 (3,3-2,z),v = 0-0 



Table 8. Line frequencies, upper energy levels and Einstein A- 
coefficients for the ortho-H20 lines which will be observed with 
Herschel/HIFI. 



Description of tine iine profiies: Most H2O lines displayed in 
Fig.[T7]have a parabolic shape, characteristic for optically thick 
unresolved emission. A few lines suffer from self-absorption in 
the blue wing (e.g., the 85^4 — 92,7 line in the ground-state). 
Particularly for lines where the optical depths at the line centre 
can be up to ~100, effective self-absorption on the blue-shifted 
side can be seen (e.g., the li — lo,i and 80,3 — 2i,2 Hnes in 
the ground-state). Lines in the first vibrational state allow one 
to trace the wind acceleration zone. These lines are considerably 
narrower than 2 f 00 (e.g., the 82,7 — 73.4 line in the j/2 — 1-state). 
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Figure 18. Comparison between H2O line profile predictions 
using (1) the Barber line list, including the ground-state and 
1^2 = 1-state (full Une, see also Fig.O, f2j the LAMDA Une hst, 
including only the ground-state (dotted line), (3) the LAMDA 
line list, including the ground-state and 1/2 = 1-state (dashed 
line), and (3) the LAMDA line list, including the ground-state, 
the 1^2 = 1- state and ly^ — 1-state (dash-dotted line). One 
can barely see the difference between including or omitting the 
1/3 = 1-state in the predictions for IK Tau, the only exception to 
this is the center of the optically thin 854 — 927 line in the ground 
state. 



Comparison to otiier line lists, including or omitting the 
1^2 and 1^3 vibrational state (see Fig. [T8]i: Using the Barber 
H2O or the LAMDA linelist yields comparable results for the 
lines displayed here. Omitting the first vibrational state of the 
bending mode (1^2 = 1), yields a decrease in line flux of a few 
percent up to 60 percent at maximum, depending on the tran- 
sition. This discrepancy will increase for lower mass-loss rate 
objects (iMaercker et alj|2008l) . The inclusion of excitation to the 
first excited vibrational state of the asymmetric stretching mode 
(1^3 ~ 1) yields a change in line flux of 2 % at maximum. This is 
in agre eme nt with the recent resu lts by iGonzalez-Alfonso et al.l 
(l2007h and IMaercker et al.l (l2009l) . who found that the inclusion 
of the 1/3 = 1-state is particularly important for low mass-loss 
rate objects. Since the Einstein A-coefficients of the symmetric 
stretching mode (i/i ~ 1) are an order of magnitude lower than 
for the asymmetric stretching, this state is not likely to affect the 
models. 



Dependence on parameters: The H2O line fluxes which wiU 
be observed with HIFI are very sensitive to certain (stellar) pa- 
rameters, and their observations will give us a handle to con- 
strain the circumstellar structure to even higher accuracy. To il- 
lustrate this, some simulations are shown in Fig. [19] 

- In case the temperature structure is approximated using a 



power law T(r) 



-0.7 



(dotted line in Fig.[T9]l, a differen- 



tial change is seen for the predicted line fluxes. Observing a 
few water lines with different excitation levels will pin down 
the temperature structure in the CSE. 

A second simulation shows the effect of using a velocity 
structure which is computed from solving the momentum 
equation (dashed-dotted line in Fig. \T% . Using another ve- 
locity law results in another gas number density and slightly 
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H2O(1,,-1„,),v = 0-0 



HjO (5,2-4,, ),v = 0-0 



HjO (3,2-22,). v = 0-0 




Figure 17. Snapshot of a few ortho-H20 lines which will be observed with Herschel/HIFI. 



different temperature structure. Using the momentum equa- 
tion, a steeper velocity gradient is obtained (see, e.g., Fig.|3]l, 
and the velocity reaches the terminal velocity at shorter dis- 
tances from the star. This results in slightly broader line pro- 
files and a flux enhancement in the blue wing since that part 
of the CSE contributing to the line profile at a certain ve- 
locity V will be shifted somewhat inward, hence attaining a 
higher source function. 

- Using a blackbody to represent the stellar radiation instead 
of a high-resolution theoretical spectrum calculated from a 
MARCS model atmosphere (see Sect. 13. Il l, only induces a 
change in the predicted line fluxes smaller than 2 % (not 
shown in Fig. \T% . For wavelengths beyond 200 fim the 
stellar flux is always represented by a blackbody in the 
GASTRoNOoM-calculations (note that the ground-state of 
ortho-water is at 23.794cm~^ or around 420 /^jm), the flux 
difference between the blackbody and the theoretical high- 
resolution spectrum around the 1^2 bending mode is shown 
Fig. [20I The reason for this negligible difference is the fact 
that the stellar radiation field is not important (in this case) 
for the H2O excitation. Excluding the stellar radiation field 
only yields a reduction of the line emission by 2 % at maxi- 
mum. 

- The dashed line in Fig. [19] shows the model predic- 
tions using the same stellar and envelope parameters as in 



iMaercker et all ( I2OO8I) : reff = 2600K, i?^ = 3.53 x lO^^cm, 
D = 300pc and M=l x IQ-^ Mg/yr. This example shows 
how maser action in the 832 — 441 transition at 621 GHz is 
very sensitive to the structural parameters. 



6. Conclusions 

In this paper, we have for the first time performed a self- 
consistent, non-LTE radiative transfer analysis on 11 different 
molecules and isotopologs C^CO, ^^CO, SiS, ^^SiO, ^OSiO, 
3"SiO, HCN, CN, CS, SO, SO2) excited in the circumstellar 
envelope around the oxygen-rich AGB star IK Tau. In contrast 
to previous studies, the temperature and velocity structure in 
the envelope are computed self-consistently, the circumstellar 
fractional abundances are linked to theoretical outer wind non- 
chemical equilibrium studies and the full line profiles are used 
as criteria to deduce the abundance structure. The Gaussian line 
profiles of HCN and SiO clearly point toward formation par- 
tially in the region where the wind has not yet reached its full 
velocity. Using the HCN line profiles as criterion, we can deduce 
that the wind acceleration is slower than deduced from classical 
theories (e.g. lGoldreich & Scovilleiri976l) . For a few molecules, 
a significantly different result is obtained compared to previ- 
ous, more simplified, studies. SiO and SiS seem to be depleted 
in the intermediate wind region due to adsorption onto dust 
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HgO (1io-1oi),v = Q-0 HgO (3q3-2i2).v = 0-Q 
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Figure 19. Comparison between H2O line profile predictions us- 
ing the Barber line list, including the ground-state and V2 = 1- 
state (1) full line: predictions using temperature and velocity 
structure as shown in Fig. |4] (see also Fig. [TTli. (2) dotted line: 
assuming a power law temperature structure T{r) ^ r^"-^, 
(3) dashed-dotted line: assuming a velocity structure consis- 
tent with the momentum equation (see discussion in Sect. 14.11 ). 
(5) dashed line: assumin g the same stellar parameters as used 
bv iMaercker et alJ (l2008h : reff = 2600K, i?^ = 3.53 x lO^^cm, 
D = 300 pc and Af = 1 x 10"^ Mg/yr. The inset in the lower right 
panel shows the full maser line profile of the 634 — 841 transition 
using these parameters. 



10'° z 




2 4 6 6 

wavelength [micron] 

Figure 20. Representation of the stellar radiation field. Upper 
panel: Comparison between a blackbody at 3000 K (gray) and a 
MARCS flux-sampled spectrum at a stellar temperature of 3000 K 
and a logarithm of the gravity of 1.5 dex (black). Lower panel: 
Comparison between a blackbody at 3000 K (gray), a MARCS 
flux-sampled spectrum (light gray) and a high-resolution theo- 
retical spectrum generated from the MARCS model with a res- 
olution of AA = 0.5A in the wavelength region between 2 and 
8 jjLva. The bandheads of the different water vibrational modes 
(the symmetric stretching mode vi, the bending mode V2, and 
the asymmetric stretching mode v^) are indicated by the vertical 
dotted lines. The inset shows a zoom around 6.3 /im where water 
absorption determines the spectral appearance. 



grains. The HCN and CS intermediate wind abundance around 
50 -300 R-i, is clearly b elow the inner wind theoretical predictions 
bv lDuari et al] ( Il999l) andlCherchneff (2006), which may either 
indicate a problem in the theoretical shock-induced modeling or 
possibly that, contrary to what is thought, HCN and CS do par- 
ticipate in the dust formation, maybe via the radical CN through 
which both molecules are formed. The lack of high signal-to- 
noise data for CN and SO prevent us from accurately determin- 
ing the circumstellar abundance stratification. It turned out to 
be impossible to model all the SO2 line profiles, particularly a 
few of the high-excitation lines. This may be due to a misiden- 
tification of the lines or to problems with the collisional rates. 
The SiO isotopic fractions point toward high ^^SiO/^'^SiO and 
^®SiO/'^"SiO ratios, which are currently not understood in the 
framework of nucleosynthesis altering the AGB isotopic frac- 
tions, but seem to reflect the chemical composition of the inter- 
stellar cloud out of which the star is born. Finally, in Sect. 15.3! 
we present H2O line profile predictions for a few lines which 
will be observed with the Herschel/HIFI instrument (launched 
on May, 14 2009). 

Appendix A: Molecular line data 

For each of the treated molecules, we briefly describe the molec- 
ular line data used in this paper Quite often, data from the 
Leiden Atomic and Molecular Database (LAMDA) are used 
dSchoier et al.ll2005l) . When appropriate, transition probabilities 
are compared to the relevant data in this database. 

CO - carbon monoxide. Both for ^^CO and ^"^CO, energy 
levels, tran sition fr equencies and Einstein A coefficients were 
taken from iGoorvit ch & Chackerian (1994). Transitions in the 
ground and first vibrational state were included up to J = 40. 
The CO-H2 collisional rate c oefficients at kinetic temperatures 
from 10 to 4000 K are from iLarsson et alJ (l2002l) . Figure [Q 
shows a good match for the transition frequencies (better than 
1 %). The Einstein A coefficients for the rotational transitions in 
the ground-state and the vibra-rotational transitions correspond 
to better than 0.5 %, but the Einstein A coefficients for the rota- 
tional transitions in the tj = 1 state may differ by up to a factor 
3 for the high-lying rotational transitions, i.e. the ones with the 
largest J quantum number. For the CO lines of interest to this 
study, i.e. rotational transitions in the w = state with Jup < 7, 
the effect on the predicted line fluxes is small: calculating the 
temperature stratification self-consistently (see Sect. 13. Il l using 
the Einstein A coefficients of the LAMDA-data base, the largest 
deviation (of 6 %) occurs for the ^^CO(l-O) Une. 

SiO - silicon m o noxid e. The SiO lineHst of 
iLanghoff & Bauschlicheil ( Il993h was used to extract the 
frequencies, energy levels and (vibra-)rotational radiative rates 
for 2^SiO, ^^SiO, and ^°SiO. Both ground and first vibrational 
state were included, with rotational quantum number J up to 40. 
The SiO-H2 collisional rates in the ground state are taken from 
the LAMDA-database. For the rotational transitions in the first 
vibrational state, it is assumed that the collisional rates are equal 
to the ones in the ground state. The vibra-rotational collisional 
rates are assumed to be zero. 

The LAMDA database only lists the frequencies and tran- 
sition probabilities for the first 40 levels in the ground state. 
Comparison with the line list of Langhoff & BauschlicheJ 
(Il993h shows that both databases nicely agree for the rotational 
transitions in the ground state (see Fig. lA.2l ). 
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Figure A.l. Comparison between the transition frequencies and 
Einstein A coefficie nts of CO as listed in the LAMP A database 
and as computed bv lGoorvitch & ChackerianI (Il994l) . 
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Figure A.2. Comparison between the transition frequencies and 
Einstein A coefficients of the rotational transitions in the ground 
state of SiO as listed in th e LAM DA database and computed by 
iLanghoff & Bauschliched (Il993h . 



SiS — silicon monosulfide. Frequencies, transition prob- 
abilities and energy levels were taken from the Cologne 
Datab ase for Molecular Spectroscopy (CDMS, iMiiller et alJ 
I2005h . Transitions in the ground and first vibrational state were 
included up to J = 40. The pure rotational transitions of silicon 
monosulfide (^^Si^^S) and its rare isotopic species have been 
observed in their ground state as well as vibrationally excited 
states by emplo ying Fourier transform microwave (FTMW) 
spectroscopy bv iMiiller et alJ (l2007h . The LAMDA database 
only lists the lowest 40 levels in the ground vibrational state, 
for which the energy levels, transition frequencies and Einstein 
A co efficients were taken from the JPL catalog dPickett et alJ 
Il998l) . The values listed in the JPL and CDMS database show 
good correspondence (see Fig. IA.3b . As for SiO, the collisional 
rates for the rotational transitions in the ground state are taken 
from LAMDA, the collisional rates in the first vibrational state 
are assumed to be equal to the ones in the ground state, while the 
vibra-rotational collisional rates are assumed to be zero. 



Rotational transitions v=0 
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Figure A.3. Comparison between the transition frequencies and 
Einstein A coefficients of the rotational transitions in the ground 
state of SiS as listed in the LAMDA database and in the CDMS 
database. 



CS - carbon monosulfide. Since the LAMDA database only 
lists the 41 lowest levels in the ground state (up to J = 40, as 
extracted from the CDMS database), the level energies, frequen- 
cies and Einstein A coefficients for the first vibrational state up to 
rotational quantum number J = 40 and the vibra-rotational tran- 
sitions between u = 1 and w = were extracted from CDMS. 
The collisional rates are taken from the LAMDA database and 
are treated as in the case of SiO and SiS. 

CN - cyanogen, cyanide radical, Transition rates 

for rotational transitions in the ground vibrational state of the 
cyanide radical are extracted from the CDMS catalog. The CN 
energy levels are indicated by three rotational quantum num- 
bers: being the total rotational quantum numbers excluding 
electron and nuclear spin, J the total rotational angular mo- 
mentum including electron spin, and F designating the spin 
quanta. All energy levels with A^up < 39 are included, yield- 
ing 235 energy levels and 508 transitions. For lack of any- 
thing better, we have applied the H2-C S collisional deexcita- 
tion rates from the LAMDA database dBlack & van Dishoeckl 
[T99I iHogeriieiide & van der Takll2000h . The typical deexcita- 
tion rates are estimated not to be in error by more than factors 
2-3 (iBlack & van Dishoecklll9"9Th . 

HCN - hydrogen cyanide. A plethora of vibrational states of 
the HCN molecule are of relevance for astronomical observa- 
tions. The states are designated by {iyii'2i^3)- The (100) mode is 
the CH stretching mode at 3311.5 cm^^; the (010) mode is the 
doubly degenerate bending mode at 712.0 cm^^; and the (001) 
mode is the CN stretching mode at 2096.8 cm^^. The excita- 
tion analysis includes radiative excitation through the stretch- 
ing mode at 3 /im and the bending mode at 14 /im. The stretch- 
ing mode at 5 /im includes transitions that are about 300 times 
weaker and is therefore not included in the analysis. In each of 
the vibrational levels, we include rotational levels up to J = 29. 
Hyperfine splitting of the rotational levels were included only in 
the J = 1 levels, where the splitting is larger than the local 
turbulent width, /-type dou bling in the 14 um tran istions was in- 
cluded. Data used are from lSchoier et al.) (l2007cl) . 
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SO - sulfur monoxide, X^Y,^. While both the JPL and 
CDMS database (and hence also the LAMDA database for 
which the values were extracted from the JPL catalog) list the ro- 
tational transitions in the v = and v = 1 vibrational state, the 
transition probabilities for the vibra-rotational transitions could 
not be found. We therefore have restricted the excitation analysis 
to the first 70 levels in the ground vibrational state as given by 
the LAMDA database. We note that the collisional rates as given 
by the LAMDA database are only listed in the range between 
50 and 350 K, and no extrapolations to higher temperatures are 
provided. In case of temperatures in the envelope being higher 
than 350 K, the rotational rates at T = 350 K were used. 



SO2 — sulfur dioxide. Energy levels, frequencies, Einstein A 
coefficients and collisional rates are taken from the LAMDA 
database (where they were extracted from the JPL database). The 
first 198 levels in the ground state are included i n the analysis . 
We note that the collisional rates are taken from lGreenl(ll995h 
and were calculated for temperatures in the range from 25 to 
125 K including energy levels up to 62 cm^^ for collisions with 
He. This set of collisional rate coefficients, multiplied by 1.4 
to represent collisions with H2, was extrapolated in the LAMDA 
database to include energy levels up to 250 cm~^ and for a range 
of temperatures from 10 to 375 K. 



H2O — water. The radiative transfer modeling includes the 45 
lowest levels in the ground state and first vibrational state (i.e. 
the bending mode 1^2 = 1 at 6.3 ^m). Level energies, frequen- 
cies and Einstein A coefficients are extracted from the high- 
accuracy computed water line list by iBarber et al.l (l2006h . cur- 
rently being the most complete water line list. E.g. while the 
LAMDA database lists 158 transition s between the 45 lo west 
levels in the ground state, the line list of lBarber et al.l (|2006|) con- 
tains 164 transitions. The Einstein A coefficients for the common 
transitions agree within ^40 % (see Fig. lA.4l l. 

The H2O-H2 collisional rates in t he ground state are taken 
from the H20-He rates bv iGreeiTe t al. (1993), corrected by 
a factor 1.348 to account for collisions with H2. Rotational 
collision rates within the first excited state are taken to be 
the same as for the ground state, while collisions between 
the ground and first excited state 1^2 are based on the 
ground state rotational collision rate c oefficie nts scaled by a fac- 
tor 0.01 (Deguchi & Nguven-Q-Rieul [T990l) . In their analysis, 
iDeguchi & Nguven-O-Rieui (il990i) show that the uncertainties 
in the vibrational colUsional rates have no effect on the calcu- 
lated H2O lines. 



^tntionol tronsitloris v = :\ 



In Sect. 15.31 we compare the H2O line profile predictions us- 
ing the Barber and LAMDA line list and the effect of excluding 
the 1^2 = 1 bending mode and including the asymmetric stretch- 
ing mode = 1 is discussed. 
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Figure A.4. Comparison between the transition frequencies and 
Einstein A coefficients of the rotational transitions in the ground 
sta te of H2O as listed i n the LAMDA database and as computed 
bv lBarberet"an (120061) . 
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